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Abstract 

With the increase of complexity and coherence of superconducting systems made using the principles of 
circuit quantum electrodynamics, more accurate methods are needed for the characterization, analysis and 
optimization of these quantum processors. Here we introduce a new method of modelling that can be 
applied to superconducting structures involving multiple Josephson junctions, high-Q superconducting cav¬ 
ities, external ports, and voltage sources. Our technique, an extension of our previous work on single-port 
structures [l| , permits the derivation of system Hamiltonians that are capable of representing every feature 
of the physical system over a wide frequency band and the computation of Ti times for qubits. We begin 
with a “black box” model of the linear and passive part of the system. Its response is given by its multiport 
impedance function (w), which can be obtained using a finite-element electormagnetics simulator. The 
ports of this black box are defined by the terminal pairs of Josephson junctions, voltage sources, and 50H 
connectors to high-frequency lines. We fit (w) to a positive-real (PR) multiport impedance matrix 
Z (s), a function of the complex Laplace variable s. We then use state-space techniques to synthesize a finite 
electric circuit admitting exactly the same impedance Z (s) across its ports; the PR property ensures the 
existence of this finite physical circuit. We compare the performance of state-space algorithms to classical 
frequency domain methods, justifying their superiority in numerical stability. The Hamiltonian of the mul¬ 
tiport model circuit is obtained by using existing lumped element circuit quantization formalisms @,0. Due 
to the presence of ideal transformers in the model circuit, these quantization methods must be extended, 
requiring the introduction of an extension of the Kirchhoff voltage and current laws. 
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1. Introduction 


Superconducting electronics is one of the most promising candidates for the realization of the hardware 
of a quantum computer. Small scale superconducting quantum mocessors have been demonstrated using 
coplanar waveguide (CPW) circuits and 3D microwave cavities [J, la,l6|. They typically consist of Josephson 
junctions coupled to cavity resonators, which in turn are coupled to microwave feedlines that provide the 
functions of readout and quantum gate operations. Josephson junctions are lossless nonlinear circuit elements 
providing the anharmonicity needed to have a quantum energy spectrum in which two unique low-lying 
energy levels can be picked out to define the qubit. 

From a classical point of view, superconducting quantum processors are microwave systems having multi¬ 
ple resonant modes. The significant increase in quality factors (Q factors) of both qubit and cavity modes in 
the last decade requires highly accurate models for the design and optimization of those systems. Previous 
approaches to model multi-mode superconducting qubit systems encountered problems in estimating loss 
rates Q- Our work has been stimulated by a recent study Q which introduced the “black box" concept 
for separating the modeling of these systems into the linear, passive, distributed part, connected via ports 
to the nonlinear or active parts of the system. We have improved upon the circuit modeling used in @>[11) 
which is based on the so-called Foster ^proach 0 which is not fully justified for lossy systems. 

In our first work in this direction (l| we have introduced a method for the accurate characterization 
of superconducting microwave circuits involving a single Josephson junction connected to a one-port black 
box. We have shown that the ad-hoc extension of Foster synthesis can be replaced by the exact impedance 
synthesis technique of Prune 0, which permits the full response of the black box, both reactive and 
lossy, to be matched to very high accuracy over a large frequency range. Applying previous formalisms ii 
developed for the quantization of lumped element circuits, we have shown that one can derive highly accurate 
Hamiltonians for single qubit systems. We have also estimated relaxation rates of the qubits, showing that 
the results are systematically improved compared with Q . 



Figure 1: Multiport Blackbox Impedance. The port electromagnetic response of a linear and passive structure is represented 
by a multiport blackbox having impedance Z(s). The N ports of this black box can be shunted by multiple Josephson 
junctions Lj^s, voltage sources V or resistors R^. In this paper we present systematic techniques for deriving a highly accurate, 
lumped-circuit description of this black box response, and we obtain the system Hamiltonian and relaxation rates for the full 
system. 
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In the present paper, we extend our method to handle the task of modeling and quantizing black-box 
systems involving multiple ports, with multiple Josephson junctions along with voltage sources and external 
impedances. We again follow the black box approach by representing the linear and passive part of the 
microwave circuit by a multiport impedance Z (s) (s is the complex Laplace variable) as shown in Fig. ([I} 
with the Josephson junctions, voltage sources and resistors shunting the ports of this multiport structure. 
The first step in the modeling is to perform a synthesis of the multiport impedance matrix. Synthesis is 
a technical term from theoretical electrical engineering, meaning the systematic determination of a finite 
lumped element electrical circuit which admits the impedance matrix Z (s) across its ports. Often, synthesis 
is to be done approximately, or within certain design specification. Here the goal is simple: the synthesis 
should be exact. With certain caveats, we do indeed achieve such an exact s^thesis. With the synthesized 
circuit, we can then apply lumped-element circuit quantization formalisms [3) 121 (with some new adaptations 
to be described in this paper) to derive a Hamiltonian for the circuit and to compute relaxation rates. 

To obtain the impedance Z (s) of a physical structure, one can imagine a variety of approaches. Ideally 
this function would be extracted from experiment, but it is hard to probe (e.g., by a spectrum analyser) the 
response at the Josephson-junction ports, which are inaccessible to external noninvasive contacting. But it 
is felt that electromagnetic response calculators are quite reliable for reliably simulating this response. The 
simulation involves the drawing of a 3D model of the structure to be represented in the blackbox in Fig. © 
using a finite-element simulator such as HFSS [l^ . In such programs, ports may be conveniently defined 
across the metal nanostructures defining the Josephson junctions, or where the coaxial inputs enter the 
cavities as shown in Fig. (USD. The simulated impedance response hsim (w) of the system is then computed 
by solving Maxwell’s equations at a discrete set of frequencies over a finite frequency band. If N ports are 
defined Z^im (w) will be a x matrix, with entries being complex numbers at every frequency. 

Next we fit 'Lsim (w) to an impedance matrix Z (s) with a finite number of poles Np. For that purpose, 
an extensively developed formalism and software package known as Vector Fitting (VF) is available, and we 
make use of it here [l3j . VF takes the number of poles Np as input and outputs an impedance matrix Z (s) 
with Np poles as in Eq. (d© which is a least-squares fit to hsim (w) over the given frequency band. The 
Rfc’s are the residues of the poles at finite (complex) frequency Sk , and E is the residue at infinity. Passivity 
of Z (s) is enforced by a subroutine of VF (l^ . 


Z(s) = 


Np 


k=l 


s — Sk 


D -f Es 


( 1 . 1 ) 


Z (s) must be Positive-Real (PR) (III . Il^ for a finite physical lumped element circuit to exist having the 
exact multiport impedance Z (s) across its ports. We have recapitulated the PR conditions in the s-domain 


for a one-port impedance function z (s) in (l|. In Section (12.41) we state the PR conditions in state space for 
a multiport impedance Z (s). Z (s) generated by VF is PR since VF produces stable poles and enforces the 
passivity. 

Below we start by revisiting the one-port problem from the state-space perspective. The state-space 
approach is central to modern circuit modeling and to control theory in a wide range of engineering disciplines 
[IGl Il7l|. In short, this approach represents the dynamics of the system in the time domain, identifying a 
sufficient set of variables so that the full dynamics is described by differential equations that are first-order 
in time. In Section (12.21) we give an introduction to the state-space theory, keeping it at a level that will 
be enough for our purposes. Working in state space allows for numerically more stable circuit synthesis 
algorithms and leads to a straightforward aRorithm for the generation of multiport circuits treatable in 
existing circuit quantization formalisms [3, Q, when extended in the way described shortly. In Section 
(12.41) we re-state the PR conditions in state space and in Section (12.51) we describe the extension of Brune’s 


method of circuit synthesis in state space m for one ports. We then show how to quantize the one-port 
state-space Brune circuit in Section (12.61) with the help of an “effective Kirchhoff” technique which eliminates 
ideal transformers from the circuit equations, leading to equations that can be thought of as generalizations 
of the Kirchhoff current and voltage laws. In the main part of this paper in Section ([3|), we extend our 
analysis to multiport circuits. We describe the multiport Brune algorithm in state space in Section (13.11) 
and apply the effective Kirchhoff method to quantize the multiport Brune circuit and compute loss rates 
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of qubits in Section (13.21) . We consider one-by-one the cases of the multiport Brune circuit being shunted 
by Josephson junctions in Section p.2l) . external resistors in Section (13.31) and voltage sources in Section 
(EH- We show the explicit application of our method on a 2-port 1-stage example circuit in Section dm 
and study a numerical example in Section (Ol) corresponding to a realistic 3D transmon setup, which is a 
3-port circuit for which the Brune synthesis gives a 12-stage circuit. 

2. One-Port Brune Circuit Quantization in State Space 

2.1. Introduction 

In Section II of [ij we presented Brune’s algorithm in s-domain (or frequency domain), s being the 
complex Laplace variable. In s-domain the impedance is given as a rational function of s (or as a rational 
matrix for multiports). The Brune algorithm in the s-domain requires a partial-fraction expansion for 
the determination of circuit parameters. Partial fraction expansion is an ill-conditioned operation since it 
requires finding roots of polynomials. Root finding is a numerically unstable problem: the roots become 
very sensitive to small perturbations in the coefficients as the degree of the polynomial increases. This is 
illustrated by Wilkinson’s polynomial [l^ . The problem becomes even more severe if one wants to apply 
multiport generalizations of synthesis algorithms (l5| since the degrees of polynomials increase with the 
number of ports. Synthesis methods given in s-domain are usually referred to as classical network synthesis. 
See [lij for a comprehensive summary of classical synthesis algorithms. Classical synthesis methods appeared 
first in the historical development of the subject and played a key role in building the theory and expressing 
synthesis procedures. However they are not suitable for computer implementation due to the stability issues 
mentioned above. 

The situation however is not so hopeless. Network synthesis algorithms can also be expressed in the state 
space. The state-space approach can be seen as a reformulation of the synthesis problem in the time-domain. 
See [I^ for a comprehensive coverage of network theory in state space. Most synthesis methods reexpressed 
in state space requires the solution of a type of Riccati equation 0- Solving the Riccati equation when 
the system’s poles approach the imaginary axis is a numerically “hard” problem. Since superconducting 
circuits have very little loss we usually encounter hard instances of Riccati equations in our models. Brune’s 
algorithm expressed in state space provides a method for solving such hard Riccati equations. Reducing the 
complexity of the problem by a small amount at each step it avoids numerical instabilities appearing in more 
direct methods which try to complete the synthesis in fewer steps. See for a discussion of how Brune’s 
method in state space might help in solving hard Riccati equations. In the following section we briefly review 
state-space formalism and present Brune’s impedance synthesis algorithm expressed in state-space terms. 

In Section (12.611 we introduce a new technique for the quantization of the one-port Brune circuit with 
ideal transformers. We call this new technique the “effective Kirchhoff” method. We will first see how much 
the effective Kirchhoff method simplifies the analysis of the Brune circuit presented in [ij. However the 
full power of the effective Kirchhoff method will be apparent in the last section of this paper when we will 
apply it to quantize the multiport Brune circuit. The earliest appearance (and the only one that we could 
find) of this technique in the literature is in where ideal transformer variables are eliminated from mesh 
equations to compute some effective mesh impedance matrices. 

2.2. State-Space Formalism 

In the state-space formalism (^e Chapter 3 of [I^ for more details on the state-space formalism in the 
context of network synthesis or [iTj in the context of dimensionality reduction theory) the state of a linear 
time-invariant system with m inputs and n outputs is given by a real vector x of length N. The time 
evolution of the state is described by a first-order differential equation 

X = Ax -I- Bu (2.1) 

where u is the input vector of length to, A a (TV x N) matrix and B a (TV x to) matrix. The output vector 
y is related to the input vector u by the following algebraic relation which involves also the state vector x 
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y = Cx + Du (2-2) 

The output vector y is of length n, C is a (u x N) matrix and D a (n x m) matrix. If u holds the currents 
and y holds voltages at the ports of a network then m = n and the multiport impedance is given by 

Z(s) = D + C(sl-A)”^B (2.3) 

We will only consider real realizations here such that the matrices {A, B, C, D} are all real. 

Now let’s assume that we transform the state x by a non-singular transformation T such that the new 
state xi is given by 


Xi = Tx 

Then using Eqs. m and (1^ the state-space description for the state xi is given by 


(2.4) 


where 


xi = AiXi + Biu 

(2.5) 

y = CiXi -I- Diu 

(2.6) 


Ai 

= TAT’i 

(2.7) 

Bi 

= TB 

(2.8) 

Cl 

= CT-i 

(2.9) 

Di 

= D 

(2.10) 


The important point to note here is that the input-output relationship is unchanged that is {Ai, Bi, Ci, Di} 
is another state-space realization for the impedance Z (s), if u and y are the currents and voltages at the ports 
of the network corresponding to Z (s), respectively. To show this let Zi (s) be the impedance corresponding 
to the realization {Ai, Bi, Ci, Di} then by Eq. (12.311 


Di-t-Ci(sI-Ai)"^Bi 

(2.11) 

D + CT-^ (si - TAT-i)”^ TB 

(2.12) 

D + C(sl-A)~^B 

(2.13) 

Z(s) 

(2.14) 


where in the second line above we used Eos. (I2.7D-(I2.1()D. For more details see Theorem ('3.3.91 of 0I- 
2.3. Minimal Realizations 

Given the impedance Z (s), a fundamental question in state-space theory is how to find a set of real 
matrices {A, B, C, D} such that 


Z(s) = D-HC(sI-A)"^B (2.15) 

is satisfied with the dimension N of the state space being minimum. In state-space theory minimal real¬ 
izations are defined in a more abstract way. The set {A, B,C,D} is called a minimal realization for the 
impedance Z (s) if [A, B] is completely controllable and [A, C] is completely observable. 

Given the time evolution equation 


X = Ax + Bu 
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(2.16) 




[A,B] is said to be completely controllable, if given the system is at state x(to) at time to, there exists a 
control u (t) defined over [toi ^i] such that the system can be brought to the zero state x (fi) = 0 at time ti 
under the driven time evolution in Eq. (l2.1fiL 
Given the state-space equations 


X = Ax + Bu 


(2.17) 


y = Cx (2.18) 

The pair [A, C] is said to be completely observable if, given the input and ouput functions u (t) and y (t) 
over an interval it is possible to determine x(to) uniquely. 

For more details on the properties of state-space realizations and the minimal realizations we refer the 
reader to Chapters 3.3 and 3.4 of [l^ . 

For a scalar impedance function z (s) the problem of finding a minimal state-space realization is relatively 
easy to answer. Without loss of generality we can assume that z (oo) = 0. Let z (s) be given as 


^ 62S -f bi 

s" -I- -I- ... -I- 023 + Oi 


(2.19) 


We assume that the numerator and the denominator polynomials in Eq. (I2.19|) have no common factors. If 
we define A in companion matrix form as 


/ 0 1 0 

0 0 1 


A = 


0 0 0 
\ -tti -02 -03 


^ \ 

0 

1 

/ 


( 2 . 20 ) 


together with the following definitions for B and C 


f ^ \ 
0 


B = 


0 

V 1 / 




/ &i \ 
^2 


^n — 1 

\ K / 


( 2 . 21 ) 


( 2 . 22 ) 


then {A, B, C} is a minimal realization for 2 : (s). This is equivalent to [A, B] being completely controllable 
and [A, C] completely observable. 

Finding a minimal realization corresponding to a multiport impedance matrix Z (s) is more involved and 
there are many procedures to find one. We will follow a physically motivated approach to find a minimal 
realization. The fitted impedance Z (s) obtained by VF in Eq. (ED contains most of the time numerical 
noise which makes its residue matrices full rank. This is generically unphysical since a full rank residue 
matrix would correspond to a degenerate mode at a finite frequency. Finding a minimal representation for 
such an impedance would introduce unphysical degrees of freedom. To cure this problem we will apply the 
“compacting” technique described in [22] to reduce the rank of residue matrices and to obtain a minimal 
realization for our models. 
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Figure 2: Brune circuit extraction in state space. 


Model-order reduction techniques are also used to reduce the dimension of non-minimal realizations. In 
applying order reduction procedures one should make sure that the passivity and reciprocity of the system 
is preserved (l^. 1^. 

2.4- Positive-Real Property in the State Space 

Given an impedance function z (s) (or an impedance matrix Z (s)) an important question is whether it 
corresponds to a finite physical circuit. For the one-port case Brune answered this question by introducing in 
[HI the so called “Positive Real (PR)” conditions. There exists a finite physical circuit having the impedance 
z (s) across its terminals if z (s) satisfies the PR conditions. Here we state PR conditions given in [l| in 
frequency domain for a one-port impedance function in the state-space language for the most general case 
of a multiport impedance matrix Z (s). 

Positive Real Lemma 

Given an mx m impedance matrix Z (s) corresponding to an m-port network with Z (c») < oo and with 
a minimal realization {A, B, C, D}. Z (s) is positive real if and only if there exist real matrices P, L and 
Wq with P being positive definite symmetric satisfying 


PA + A^P 

= -LL^ 

(2.23) 

PB 

= - LWo 

(2.24) 

WjWo 

= D + D'^ 

(2.25) 


Popov Lemma” 

11,0,113,Hi. For more 


The Positive Real Lemma stated above goes also under the name “Kalman - Yakubovich 
in control theory literature which refers to names involved in its development [1 
details on Positive Real lemma see Chapter 5 of [l^ . 

Most of the synthesis algorithms stated in state space Q are based on the determination of the matrix 
P which is usually done by solving a Riccati equation. The algorithm we present in the following will ide ntify 
P in a recursive way which avoids numerical difficulties appearing in more direct methods presented in |l6l |. 

We will now present the Brune algorithm in state-space terms as described in [I3| . 


2.5. Brune’s Algorithm in the State Space (One-Port Case) 

Here we assume that we have a one-port positive real impedance function z (s) with the minimal real¬ 
ization {A, B, C, D} (We note that D = D is a scalar in this case). As shown in Fig. Q Brune’s algorithm 
in state space starts with the extraction of the series resistance i?i. Using Eq. (12.31) we can evaluate the 
real part of the impedance over the imaginary axis as follows 
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Re[z (jw)] 


for some frequency ojq with 


i {z {juj) + z i-juj)) 

(2.26) 

D+^C {jojI - A)-i B + ic {-jujl - A)-i B 

(2.27) 

D-CA (w^I-f A2)"^B 

(2.28) 

given by 


Ri = min Re[z (jui)] 

UJ 

(2.29) 

Re[z (jwo)] = i?i 

(2.30) 


Let the network N 2 in Fig. m be described by the state-space equations 


X2 = A2X2-t-B2U2 (2-31) 

y2 = C 2 X 2 + D 2 U 2 (2.32) 

so that the realization {A 2 , B 2 , C 2 , 1 ? 2 } corresponds to the impedance Z 2 (s) = £>2 + C 2 (si — A 2 )~^B 2 
seen at the terminals of the network N 2 {D 2 is a scalar). 

Then the state-space equations for the network Ni are given by 



(2.33) 


where xci = 3'iid xli = ■ Hence the state-space equations for the network Ni are of the 

form 


with 


xi = AiXi-l-BiUi 

yi = CiXi-kDiUi 


Xi = ( x|’ XCi XLi )' 


Ai = 


/ A 2 
0 




Til V LiCi 

C? _ 1 _ D 7 . 


Bi = 



(2.34) 

(2.35) 

(2.36) 

(2.37) 

(2.38) 


























_ f C 2 (1 1/ni) 

’-"1 - I, m Vc: 


D 2 

nlVLi 


(2.39) 


D, = 


Dn 


(2.40) 

Ui 

The realization {Ai, Bi, Ci, Di} then corresponds to the impedance function zi (s) = Di+Ci (si — Ai)~ Bi 
seen at the terminals of the network Ni (Di is a scalar) which is related to z (s) by 

zi (s) = z (s) - i?i (2.41) 

We note that Re [zi (jwo)] = 0. 

The following lemma stated in [l^ shows that if zi (jcuo) + zi (—JWq) = 0 is satisfied for some wq > 0 
for a positive-real impedance function Zi (s) with a minimal realization {Aq, B^, C^, Da} then there exists 
a coordinate transformation T which would give an equivalent state-state description {Ai, Bi, Ci, Di} for 
the impedance Zi (s) in the form given in Eqs. (12.34112.4111) with 


Ai = TA,T 
Bi = TB„ 

Cl = CaT-^ 
Dl = Da 


-1 


(2.42) 

(2.43) 

(2.44) 

(2.45) 

See Section (12.211 for why {Ai,Bi,Ci , Dl} is an equivalent realization for the same impedance zi (s). 

An explicit procedure is presented in to compute T. We now state the lemma and describe the 
algorithm to compute T. 

The Fundamental Lemma (one-port case) 

Let zi (s) be a positive-real impedance function with the minimal realization {Aq, B^, C^, Da} satisfying 
zi (jwo) + zi (—jwq) = 0 for some finite frequency wq (with jwo not being an eigenvalue of A^). Then there 
exists a coordinate transformation matrix T such that Ai = TA^T'^, Bi = TB^, Ci = and 

Dl = Da are of the form given in Eqs. (12.34112.4lJl) . 

Now we show how to construct T. 

1) Construct a nonsingular matrix with the last two columns of being (wgl + Ba and 

-A„ (w2i + a2)-'b„. 

2) Set Ab = ToAoT^ \ Bf, = T^Ba and C;, = ^ and compute 


Cf, (ujH + AI 


-1 


Cb (wgl + Al) Ab 

where K 22 is a 2 x 2 matrix. Define 


= K 


12 


K 


22 


(2.46) 


' ^22^12 I 

3) Set Ac = TbAfjTjj Be = T^B^ and Cc = CbT^b ^ fhon 


for non-zero a, fd. Define 


ThenT = TcT6Tc. 


Cc Kl + A2j 

:cKi + a2) ^a 


-1 

-1 


Tc = 



0 0 
0 0 ^2 


(2.47) 


(2.48) 


(2.49) 
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Figure 3; The capacitive degenerate Brune stage corresponding to the extraction of the resistor Ri in Eq. 112.291 1 at inifinite 
frequency tug = oo. 


2.5.1. The Capacitive Degenerate Stage 

It is possible that the frequency ujq in Eq. (12.301) where the minimum in Eq. (12.291) is reached occurs at 
infinity ujq = oo. In such a case we need to extract a capacitive degenerate stage which doesn’t involve the 
inductive circuit as shown in Fig. ([2|). (The case wq = 0 requires the extraction of an inductive degenerate 
stage as shown in Fig. 0, see Section (12.5.21) below). 

Let the network N 2 in Fig. m be described again by the state-space equations 


X2 = A 2 X 2 -1- B 2 U 2 (2.50) 

y2 = C 2 X 2 -I- D 2 U 2 (2-51) 


for some real matrices {A 2 , B 2 , C 2 , D 2 } so that the realization {A 2 , B 2 , C 2 , II 2 } corresponds to the impedance 
Z2is)=D2 + C2{sl-A2)-^B2 seen at the terminals of the network N 2 . 

Then the state-space equations for the network Ni are given by 



yi 



(2.52) 

(2.53) 


where xci = and D 2 is a scalar in the one-port case. Hence the state-space equations for the 

network Ni are of the form 


with 


Xi = AiXi-l-BiUi 

(2.54) 

yi = CiXi+DiUi 

(2.55) 

Xi = ( x|’ XCi 

(2.56) 

/ A B 2 C 2 B 2 \ 

= ’ c, '’■‘F 

(2.57) 

\ D 2 y/Ci D2C1 J 


w 

II 

(2.58) 


II 




























Figure 4: The inductive degenerate Brune stage corresponding to the extraction of the resistor Ri in Eq. 112.291 1 at zero 
frequency u;o = 0. 


Cl = ( 0 ^ ) (2.59) 

£>1 = 0 (2.60) 

The realization {Ai, Bi, Ci, £> 1 } then corresponds to the impedance function zi (s) seen at the terminals 
of the network Ni which is related to z (s) by 


Zi (s) = z (s) - i?i 


(2.61) 


In such a degenerate case we should also modify the Fundamental Lemma as follows: 

The Fundamental Lemma (one-port capacitive degenerate case) 

Let zi (s) be a positive-real impedance function with the minimal realization {A^, Bq, Cq, Da} satisfying 
Zi (juJo) Zi (—jujo) = 0 for ujQ = 00. Then there exists a coordinate transformation matrix T such that 
Ai = TAaT“^, Bi = TBa, Ci = and £>i = Da = 0 are of the form given in Eqs. 

Now we show how to construct T. 

1) Construct a nonsingular matrix with the last column of being B^. 

2) Set Ab = TaAoT^^, Bf, = T^Bq and Cb = and make the partitioning 


Cb = ( Ki 2 K22 ) 

where K 22 is a scalar. Define 

( K22IK12 I ) 

3) Set Ac = ThAftTji Be = T^Bf, and Cc = C^Tj^ ^ then 

Cc = ( 0 ••• 0 ) 

for a non-zero a. Define 


Then T = TcTfeT^. 


Tc 


I 

0 


0 

a 


(2.62) 


(2.63) 


(2.64) 


(2.65) 
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Figure 5: Brune circuit obtained from the application of the state-space Brune algorithm is shown in the dotted box. We note 
that the coupled inductors at each stage of the classical Brune circuit in Fig. (1) of Q are replaced by ordinary inductors 
shunting ideal transformers as shown in Fig. I©. 


2.5.2. The Inductive Degenerate Stage 

It is possible that the frequency wq in Eq. (12.301) where the minimum in Eq. (|2.29l) is reached occurs at 
zero, Wo = 0. In that case one needs to extract an inductive degenerate Brune stage as shown in Fig. 

It is straightforward to modify the state-space Brune algorithm to synthesize an inductive degenerate stage 
in a way similar to the treatment we did for the capacitive degenerate case above; we do not present this 
algorithm here. See Section (12.911 for how one can treat such a degenerate stage in the circuit quantization 
and dissipation analysis with the assumption that the loss introduced by the resistor Ri in Fig. (U) is small. 

2.6. Quantization of the One-Port State-Space Brune Circuit 

We call “the state-space Brune circuit” the circuit obtained by the application of the state-space Brune 
algorithm described in the previous section. An M stage state-space Brune circuit is shown in the dotted 
box in Fig. We note that the coupled inductors at each stage of the classical Brune circuit in Fig. (I) 
of [ij are replaced by ordinary inductors shunting ideal transformers in Fig. which is justified by the 
equivalence shown in Fig. ®. 



Figure 6: Equivalence of the inductive circuit generated by the state-space Brune algorithm to the tightly-coupled 
pairs appearing at each stage of the original Brune circuit in Q. Circuit parameters are related by L = L\\ and n = 

Ideal transformers were not in the toolbox of any previous circuit-quantization analysis i,i 
treat them we will introduce here a new technique which will eliminate them by generating effective loop 
matrices involving turns ratios in their entries. We will see how this technique will simplify significantly the 
analysis of the one-port Brune circuit presented in [ij. It will allow us to skip the transformation dehned in 
Eq. (A9) of [l|. We will however see the full power of this technique in Section (13.21) when we will use it to 
quantize the multiport Brune circuit. To analyze this circuit we need to modify it to make it treatable in 
the formalism of |3|. 


inductor 


to 
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Figure 7: Modified state-space Brune circuit. Tree branches are shown in black and chord branches are shown in blue. Formal 
capacitance Cm+i introduced for a technical reason: with the substitution Cm+i = able to compute 

dissipation rate due to shunt resistor Rm+i in the formalism of Ql- After the coordinate transformation (see below) we take 
Cj —)■ 0 limit. 


An augmented form of the state-space Brune circuit is shown in Fig. ([7]). The last resistor Rm+i is 
replaced with a capacitor Cm+i which is included in our analysis later through the substitution Cm+i 
1/{iuiRm+i)- Its contribution to the dissipation rate will be computed referring to the equation of motion 
Eq. (61) in 0. 


1 : ni 



1 : um 



Cm 


Figure 8: Lossless part of the state-space Brune circuit. It is this circuit that corresponds to the system Hamiltonian derived 
below. As discussed above, we take the limit Cj —>■ 0 so that this element is removed. The lossless circuit is obtained from 
Fig. 0 by taking Ri, ■■■Rm 0 and Rm+i It is these different limiting treatments that require the descriptions of 

the in-series resistors Ri — Rm follow the low-impedance treatment as in @|, while the description of the shunt resistor Rm+i 
needs the high-impedance treatment as in [2|. 

The lossless part of the state-space Brune circuit which corresponds to the system Hamiltonian derived 
below is shown in Fig. ([5]). As shown in Fig. (O in the special case of unity turns ratio, this circuit is one 
of the lossless Foster forms in [l^l- We again add a formal capacitance Cj shunting the Josephson junction. 
This is required for a non-singular capacitance matrix if there are no degenerate stages. 


1 : 1 



Figure 9: Circuit identity showing that inductor-ideal transformer pairs appearing in state-space Brune stages simplify in the 
case of turns ratio equal to one; in this case Fig. becomes identical to one of the classic lossless Foster canonical forms in 

na. 


We now show how to treat ideal transformers by extending the loop analysis in 0 . Kirchhoff’s laws are 
given by Eqs. (4-5) in 0 


FIch = —Itr 


( 2 . 66 ) 
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F^Vtr = Veh 


(2.67) 

where we have assumed that there are no external fluxes in circuit loops. F is the loop matrix with entries 
being 0, 1 or —1 derived by a graph theoretical analysis of the circuit Q. After the effective Kirchhoff 
analysis done below F will be replaced by the effective loop matrix with real-valued entries. Itr and 
Ich are the tree and chord branch current vectors respectively partitioned as follows 

Itr= (2.68) 

Ieh= (2.69) 

Here labels J, L, Z, C, T correspond to Josephson junction, inductor, resistor, capacitor and ideal trans¬ 
former branches, respectively. and are the current vectors for the ideal transformer branches in 
the tree and chords respectively. We also partition loop matrix F according to the partitioning of current 
vectors 


Fjc 

F JT \ 

F LC 

F lt 

F zc 

F ZT 

Ftc 

Ftt } 


(2.70) 


We will eliminate ideal transformer branches from Kirchhoff laws in Eqs. (I2.66l) - (|2.67l) to get an effective 
loop matrices and (F^)^ ^ such that we have a new set of effective Kirchhoff relations 


■peZ/je// 
^ -^ch 


re// 


(FT)«//V«// ^ ^e// 


(2.71) 

(2.72) 


where 


I 


eff 

tr 


[IjM.Iz) 


I 


eff 

ch 


ic 


and 


F®// = 


(jTyff ^ 



(2.73) 

(2.74) 


We note here that the entries of the effective loop matrix are real numbers, being functions of ideal 
transformer turn ratios as we will see below. 

In this section for simplicity reasons we will derive only the effective Kirchhoff’s current law in Eq. (12.711) 
by computing . We postpone the derivation of the effective Kirchhoff’s voltage law and the computation 
of the matrix to the Appendix (17.2.11) . However we note here that 

(jTyff ^ (jefyT 75 ) 

should be verified to hold. This ensures the symmetry of various matrices computed in the formalisms of 
like the capacitance matrix C and the stiffness matrix Mq for example. 
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Now we claim that in Eq. (12.73^ is given by 


n{/ = ( 1 1 11) 


(2.76) 


/l (1-ni) ••• (1-ni) i^-ni) \ 


■pe// _ 
^ LC — 


VO 


1 (l-riM-i) (l-riM-i) 

1 (1 - um) J 


(2.77) 


pe// _ 

^ zc — 


/ 1 1 
1 


1 

1 1 


VO 11/ 

where is a row vector of length (M +1), and are M x (M + 1) matrices. 


matrices in Eqs. (I2.76I1 - (I2.78I) can be computed we first note the following 

iP = -FtcIc 


with 


Ftc = 


/ 0 1 1 
0 1 


1 \ 


Vo Of/ 

where Ftc is a M x (M + 1) matrix. We note that Ftc doesn’t involve any turns ratios. 

(ch) (tr) 

transformer relations = — NI^ with N being the diagonal matrix of turns ratios 


/ 


N = 


ni 


V 0 


Um 


and Eq. (I2.79|) we get 

Inductor currents are given by 

where 


= NFtcIc 


1l = -FlcIc - 


F LC = 


(1 1 
1 


V 0 


1 1 
1 1 


\ 

1 1 ) 


(2.78) 

To see how the 

(2.79) 

(2.80) 

Using the ideal 

(2.81) 

(2.82) 

(2.83) 

(2.84) 


and Ftt = —I- Using Eqs. (I2.82|l and (12.831) we get 


II = — {F LC — NFtc) Ic 

which gives the effective loop matrix F^/^ 


(2.85) 
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( 2 . 86 ) 


Flc = Flc - NFtc 

/I (1-ni) • 


VO 


(l-ni) \ 

{l — UM-i) (1 —riM-i) 

1 (I-tim) / 


(2.87) 


We note 


F 


e// 

JC 


that F^-^ is no longer 
is simply given by 


a binary matrix as 


we have turns 


ratios appearing in its entries. 


( 2 . 88 ) 

(2.89) 

Since the current through the Josephson junction depends only on chord capacitor currents 


Fyj = Fjc 

=(11... 11) 


Ij = -Fjclc (2.90) 

Note that Fj{/ does not depend on turns ratios. Similarly the currents through the resistors Rj for 1 < j < 
M depend only on chord capacitor currents 


Iz = —Fzc^c 


Hence 


re// _ 

zc ~ 


F zc 

( 1 1 

1 

V 0 


(2.91) 


1 M 
1 1 


1 1 / 


Now one can write an equation of motion for the one-port state-space Brune circuit in Fig. dSD in the 
form of Eq. (29) in Q 


(C+Cz)*i = - 


9$ 


(2.92) 


However effective loop matrices derived above have to replace the ordinary loop matrices while computing 
quantities appearing in the equation of motion Eq. (|2.92ll like the capacitance matrix C and the dissipation 
matrix C^, as we now show. 

We compute the capacitance matrix Cq for the Brune circuit in Fig. ([S]) using the Eq. (22) of Q with 
the effective loop matrix 


Cn = 


Cj 0 

0 0 




where C is the diagonal matrix of capacitances 

/ 


c = 


(2.93) 


(2.94) 


V 0 


Cm+1 


17 









and 


■peff _ 

•^c — 


JC 
-pe// 
^ LC 


Lt in 


Eq. (15) of 0 is a diagonal matrix of inductances 


(2.95) 


With 



(2.96) 


(2.97) 


since there are no 


chord inductors. We get using Eq. (31) of 0 


Mo = 

'0 0 

.0 Lr 


(2.98) 

(2.99) 


We skip the first transformation defined in Eq. (A9) of [ij and the truncation afterwards since we 
are already in the low dimensional subspace with (M + 1) degrees of freedom. We again define a local 
transformation matrix T which makes the Langrangian description (i.e., both Co and Mo) of the system 
band-diagonal: 


/ 1 

-1/(1-ni) -1/(1-ni) 

1 /( 1 -7^2) 1/(1-712) 


0 




( 2 . 100 ) 


V 0 

Applying T to Co and Mo we get 


(-l)^/(l-„^) /(i_nM)) 


C = T*CoT 

(Cj + nlC[ mC; 
niC^ Cl + n2C2 

0 


0 


^M-1 “C 


UmC 


M 


nMCju + Cm+i/ 


( 2 . 101 ) 


( 2 . 102 ) 
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(2.103) 


Mo = r‘MoT 



where C'- = Q/ (1 - n^f, L'^ = Lj{l- njf. 

A Lagrangian Cq (and equivalently a Hamiltonian 'Hs) can be written as 

Co = -u{^),ns = \q^c-^Q + u (^) 

where 


(2.104) 


(2.105) 


C/(^) = -(^) Lj\os{^j) + ^^'^Mo^ (2.106) 

$ is the vector of transformed coordinates of length (M + 1) and d*! = (|^) <fj- We note that the trans¬ 
formation T in Eq. (12.1001) introduces a local relationship between the final coordinates $ and the branch 
fluxes ^L^’s across the ordinary inductors Lj’s for 1 < j < M in the Brune circuit in Fig. ([S]) in the 
sense that the flux is only a superposition of two consecutive coordinates and by the relation 
which gives 


for 1 < j < M. is the vector holding the fluxes across the inductors in the Brune circuit in Fig. (O such 
that 


= ^L^ ■■■ (2.108) 

with 4>,/ = (f^) fj being the flux across the Josephson junction. 

2. 7. Dissipation Analysis 

In this section our aim is to compute relaxation rates. We follow here the treatment given in [3| with 
the exception that we are going to use effective loop matrices derived in the previous section instead of the 
ordinary ones. 

We treat resistors in the Caldeira-Leggett formalism with each resistor representing a bath of harmonic 
oscillators with a smooth frequency spectrum. Couplings of the baths to the circuit degrees of freedom are 
given by ni matrices as defined in Eqs. (65) and (27) of Q and respectively. 

We will interpret the equation of motion in Eq. (29) of S as an equation of motion in Eq. (61) of [2|. 
We start by rearranging the equation of motion Eq. (29) of [3| 

= (2.109) 

Cz is given in frequency domain in Eq. (26) of Q as 
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with 


(2.110) 


Cz{oj) = mCz (w) 


m = (f%^SY 


( 2 . 111 ) 


Cz (uj) = —iujZ (uj) 



( 2 . 112 ) 


where we used Eqs. (27) and (28) of Q (correcting a typo in the sign of Cz (w)) with effective loop matrices 
and computed in the previous section. We note that is independent of ideal transformer 

turns ratios. 

Comparing Eq. H2.109I) to the equation of motion Eq. (61) in Q we identifiy in frequency domain 


Md (w) = -oj'^Cz (w) 


(2.113) 


Using Eq. (12.1101) and (12.1131) 


Md (w) = —uj'^mCz (uj) 


(2.114) 


Comparing Eq. (|2.114l) to the Eq. (64) in Q we make the identification 

L-i (w) =(w) (2.115) 

with the ni matrix being given by Eq. (|2.111l) 

= (2.116) 

where we have also taken into account the coordinate transformation T defined in Eq. (12.1001) . 

We will treat resistors one at a time [13 ■ Hence for the series resistors Rj with 1 < j < M, Cz {uj) 
defined in Eq. (12.1121) is a scalar function Czj (w) which allows us to write the kernel defined in Eq. (73) 
of [3 using Eq. (12.1151) above 


Kj (w) = (w) (2.117) 

for 1 < j < M. Also we will use the columns ihj for 1 < j < M of the matrix in defined in Eq. (12.1161) : 
nij giving the coupling of the system degrees of freedom to the bath of the resistor Rj. Eor the last resistor 
Rm+1 "we will read off the kernel Km+i (w) and the coupling vector ifiM+i directly from the equation of 
motion after the replacement Cm+i t— 1/{iuiRm+i) as shown below. 

We compute the spectral density of the bath corresponding to the resistor Rj for all resistors 1 < j < 
M + 1 using the Eq. (93) of 0 


Jj (w) = Im [Kj (w)] 


(2.118) 


where we corrected a sign typo and dropped the factor p which will be justified down below. 

Applying Eq. (124) of [3 we get the contribution to the relaxation rate from the resistor Rj {l<j<M 


1 ): 


7^ = t ■ ^1 Jj (woOcoth 
n 


f fx-um \ 

\2kBT) 


(2.119) 
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|0,1) are the qubit eigenstates of the system Hamiltonian in Eq. (I2.105P and wqi is the transition frequency 
between them. Calculating these quantities requires solving the Schrodinger equation for the system Hamil¬ 
tonian in Eq. (|2.1()5|1 above; this can be a difficult task, but many effective accurate methods have been 
developed for doing this, in many works right up to the present The vector nij represents 

the coupling of the system to the bath of the resistor Rj. 

To treat the in-series resistor Rj for 1 < j < M, we imagine that all other in-series resistors are short 
circuited such that = 0 for 1 < fc < M and k ^ j and the last shunt resistor is open circuited such that 
Rm +1 = oo. Using Eqs. (12.1161) and (12.1121) we write 

m = (2.120) 


Cz., (a;) = -zu;R, (^I + tu;R, 

where we made the replacements Z (w) ^ Rj and ^ ^'r^c being the row of F'^Jq . 

( 0 \ 


0 

entry 

(-ip+Uj+iCj+i (-lycj 

(l-nj+i) (1-nj) 


(-1)^nuCM ( —1)^ ^Cm-1 

(l-niu) ' (l-riM-l) 

\ (I-km) 


where tiij are vectors of length (M -|- 1) and 


We then have using Eq. (12.1171) 


Hence we obtain using Eq. (I2.118|) 


Cz., (u;) = 


—iujRi 


M 


1 -|- iujRj I ^ Ck 




Kj (uj) = -uj'^Czj (w) 


iuj^Rj 


( M ^ 

1 -|- iuiRj j ^ Ck 

\k=j , 


( 2 . 121 ) 


( 2 . 122 ) 


(2.123) 

(2.124) 


Jj (w) = Im [Kj (w)] 

uj^R.i 


M 


l+a;2i?2 J2Ck 


yk=j 


(2.125) 

(2.126) 
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Note that our use of the non-normalized coupling vector m_, and the flux vector $ in Eq. (12.1191) implies 
removal of the factor p (f^)^ from the dehnition of the spectral function of the bath J (w) in Eq. (93) of 
[2 (See also (|2.128|) below). 

To treat the last resistor Rm+i we replace Cm+i in the last row of capacitance matrix by 1/(iuiRM+i)- 
We get the following dissipation matrix for resistor Rm+i 


Md = Km+1 (w) mM+im^_^i, 


(2.127) 


where Km+i (w) = and tiiM+i = 


0 

VJ 


is a vector with (M + 1) rows. We then have 


Jm +1 (w) = Im [Km +1 (w)] = 




Rm+1 ’ 


(2.128) 


Compare niM+i which is dimensionless to rhj in Eq. (12.1211) which has the dimension of capacitance(Farad) 
for 1 < j < M. Note also the difference in dimensionality between Jm+i (w) in Eq. (12.1281) and Jj {+>) 
defined in Eq. (12.1261) for 1 < j < M. 


2.8. Capacitive Degenerate Case 

Here we consider only a single capacitive degenerate stage. We consider a degenerate case appearing at 
stage as described in Section (12.5.11) . In case of such degeneracy we remove the (fc + 1)*^ row in the 
matrix in Eq. (12.951) : 


/I 1 

1 (1-ni) 


tte// _ 
■^C — 


V 


0 


1 


1 


(1-ni) (1-ni) 

1 (l-rifc-i) (l-nfc_i) 

1 (1 - rifc+i) 


We also modify the T transformation defined in Eq. (|2.100p above by dropping its {k +1)*^ row to obtain 


1 1 \ 

(l-ni) (l-^^l) 

(l-nfc_i) (l-nfc_i) 
(1-nfc+i) (l-rifc+i) 


1 (1-nM-i) (l-riM-i) 

1 (I-tt-m) j 

th 


(2.129) 


/ 1 


T = 


(l-m) 


(l-n,) 

(1-112) (l-" 2 ) 




(-1)'= 


0 


(l-nfc_i) 

(l-rifc+i) ( 1 -rifc+i) 




0 


V 


(I-iim) (1-iim) / 


(2.130) 


Using again Eqs. (12.931) and (12.991) respectively we obtain after applying the modified T transformation 
in Eq. (|2.130l) above 
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0 


(Cj + n\C\ niC[ 

Ti\C-^ C-^ 112(^2 


c = 


0 


CLr 


Ck 

''^fe+iC'fc+i 


•‘k+l^k+1 




’^fe+iC'fe+i 


'•fc+2^fc+2 


'M-1 nMCf^ 

^mCj^ Cj^ + Cm+1/ 


(2.131) 


Mn = 


1 1 

L[ L[ 

A + — 

L'l L[ ^ L'^ 


1 

IX 


V 


1 

L' 


1 

IZ 


1 

lt 


0 




0 




(2.132) 


Note that the matrices above are of size M x M. 

In case of degeneracy in vectors are computed again using the Eq. (12.1201) for 1 < j < M and applying 
the T transformation in Eq. (|2.130l) . We define some auxiliary vectors 


/ 


niCi 




(- 1 ) 

(- 1 ) 


k-l rik-iCk-i 
(l-rik-i) 
k nk-^-iCk-^~l 
(1-rifc+i) 


V 

/ 


( _1 'iM—1 _2 tj\£Cm_ 

1 (1-«m) 


nif, 


0 

Cl 

(l-m) 


(- 1 ) 

(- 1 ) 


k — l Ck-i 

(l-rik-i) 
k C k+i 
(l-rifc+i) 


I Cm 

(1-nM)/ 


and the vectors nia,, and m;, 1 as 


(2.133) 


(2.134) 
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Figure 10: The inductive degenerate case. The degenerate inductor Lj is in the tree whereas the resistor Rj is replaced by the 
formal capacitor Cj which is a chord branch. Rj can be treated in the same way as the last shunt inductor Rm+i (by making 
the substitution Cj ■<— to do a dissipation analysis) in the Brune circuit provided that the loss introduced by Rj is small; 

that is a high impedance treatment is possible for Rj. 


maj (i) 


0 1 < * < j 

ma (i) j <i< M 


(2.135) 


mbj (*) 


0 1 < j < j 

nib (*) j < i < M 


(2.136) 


where nia (i), m;, (i), m.a,j (*)> ^b,j (i) are the entries of the vectors nia, nif,, m.a,j, respectively. 
Finally we define the vector liiCfc as 


(i) 


0 i ^ k 

Ck i = k 


(2.137) 


where mcj, (*) is the entry of the vector mcj,. We note that the vectors nia, nib, iiiaj, ihbj and rhc^ 
are all of length M. 

Now we can write coupling vector rhj to the bath of the resistor Rj for 1 < j < M as a function of the 
vectors defined in Eqs. (I2.135L (I2.136L (12.1371) above as 


nij =ma(j) + mfc(j) + mcfc, if j < k (2.138) 

= liia (j) + riih (j), if j > k (2.139) 

Spectral densities Jj (w) are the same as in the non-degenerate case (Eqs. (12. 12611 . (|2. 1281) 1 for all resistors. 
Note also that dissipation treatment for the last resistor Rm+i is unaffected since Cm+i is untouched in 

Eq. snm . 

2.9. Inductive Degenerate Case 

In Section (12.5.21) we noted that an inductive degenerate stage might appear in the stage of the Brune 
circuit as shown on the left part of Fig. (H®. One can deal with such a degeneracy by first replacing the 
resistor Rj extracted right before the degenerate inductor Lj with a formal capacitance Cj and putting the 
capacitance Cj in a chord branch (as shown on the right in Fig. dTUll l and treat it in the same way as the 
last shunt resistor Rm+i', that is by making the substitution Cj •<— and doing a dissipation analysis 
similar to the one done above for the last shunt resistor provided that a high impedance treatment is 

possible for Rj. The degenerate inductor Lj should be in the tree like the ordinary inductors appearing in 
the Brune circuit. We note that the degenerate inductor Lj will introduce a new degree of freedom since it 
is a tree branch compared to the degenerate capacitor treated in the previous section which didn’t introduce 
any new degree of freedom. 
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Figure 11: The multiport Brune Circuit 









































































Port 1 —!> 


Ports 


^ Ti 


n 


UlA 


+ 

yiA 


1 : 71,1 


U2A 


Li p 


l:z^i 


•ai 


+ 

y2A 


Cl vci 


UiB 

Ol bi U2_B 

+ 

; + 

yiB 

y2B 

— 

; — 


No 




N 


Figure 12: Multiport Brune circuit extraction step. The extraction starts with the Belevitch transformer Ti. The circuit that 
follows Ti is almost identical to the one-port Brune stage in Fig. lEJ except the multiport transformer coupling the first 
port to the remaining ones. 


3. Multiport Brune Quantization 

Here we extend our analysis in the previous sections to multiport circuits. The multiport Brune circuit 
is shown in Fig. m- This circuit is obtained by applying the multiport Brune’s method [l^ in state space 
which we describe below. We again apply circuit quantization formalisms ii with the help of the effective 
Kirchhoff technique we introduced for the one-port Brune circuit to derive a Hamiltonian and compute 
relaxation rates for the multiport Brune circuit. 

The circuit in Fig. CB consists of N ports and M stages. On the far left we have N terminal pairs 
corresponding to the ports. Each stage starts with the extraction of a Belevitch transformer T^. Each is 
a 2Wport transformer with N ports on the left and N ports on the right which we describe in detail in the 
next section. The circuit representation of the Belevitch transformer is shown in Fig. m- The multiport 
Belevitch tranformer is followed by a resistor extracted only at the first port. We will describe below 
how to extract Belevitch transformers and resistors. After the resistor extraction we have the reactive part 
of the multiport Brune stage. We observe that the part of each Brune stage that comes after the Belevitch 
transformer at the first port is almost identical to the one-port Brune stage. We see however an additional 
transformer izfc coupling the reactive circuit in the first port to the remaining ports, is a multiport 
transformer with the primary winding connected in parallel across the terminals of the capacitor Ck at stage 
k. Each of the secondary windings are connected in series between terminals of the Belevitch transformers 
at the remaining ports. The last stage consists of the Belevitch transformer Tm+i shunted by the resistors 
i?i,..., • 

3.1. Multiport Brune Algorithm 

In this section we will describe multiport generalization of the Brune’s synthesis algorithm described in 
state space in Section (1^ for one-port networks. Fig. (fT^ illustrates extraction of a multiport Brune 
stage. Since the algorithm is recursive we will describe it only on the first stage. At each stage the degree 
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Figure 13: The Belevitch Transformer with N ports on the left and M ports on the right. On the right side we see the detailed 

/ fll ■■■ fiiv \ 


circuit representation of this device. There is a turns ratio matrix T = 


associated with the Belevitch 


V iMl ■ ■ ■ ^MN ) _ 

transformer which relates currents and voltages on both side of the device as given in Eqs. 113.6ll and 113.711 . The arrow in the 
box is used to refer to the asymmetrical character of the Belevitch transformer. 


of the network is reduced by two hence the algorithm terminates once a constant multiport impedance is 
reached as in the one-port case. 

We will only focus on the reciprocal response case, that is when Z = Z^. We will show later how the 
gyrators appear in the multiport Brune circuit in case of a non-reciprocal impedance. 

The synthesis of a multiport Brune stage starts with the extraction of the multiport Belevitch transformer 
Ti together with the resistor ri at the first port. 

3.1.1. The Belevitch Transformer 

A generic multiport Belevitch transformer with N ports on the left and M ports on the right is shown 
in Fig. m on the left. The detailed circuit representation of the Belevitch transformer T is shown on the 
right of Fig. m which defines a M x N matrix for T 

til ■ ■ ■ tiN 
tMi ■ ■ ■ Imn 

Let the current vectors and be the vectors holding the currents at the ports on the left side and 
the right side of the Belevitch transformer T , respectively, i.e. 

(3.2) 

(3.3) 

and let the vectors and be the vectors holding the voltages at the ports on the left side and the 
right side of the Belevitch transformer T , respectively, i.e. 
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Pmn 


Figure 14: The Belevitch transformer in reflected form as used in the multiport Brune circuit in Fig. IHIL The current and 
voltage relations for this transformer are given in Eqs. lITSt . llMt . 


v(-^) = 

{vP. 

■ ■ ■, Vn J 

(3.4) 

v(-^) = 

(q*‘, 

• ■ ■ I *^M J 

(3.5) 

then we can write the Belevitch transformer relations as 



j(^J) 

— _ 


(3.6) 


v(i) = (3.7) 

One should recognize the asymmetrical character of the Belevitch transformer which we noted by putting 
an arrow in the box representing the Belevitch transformer in Fig. m- However in the multiport Brune 
circuit in Fig. (na we use the Belevitch transformer in Fig. m in a reflected form as shown in Fig. da- 

The current and voltage relations for the reflected Belevitch transformer in Fig. m are given by 

l(i) = (3.8) 

yiR) = T^yiL) (3.9) 

It is interesting to note the similarity of the Belevitch transformer relations in Eqs. dUl) and (13.71) to 
the Kirchhoff’s laws given in terms of the loop matrix F in Eqs. (12.661) and (|2.67l) given in Section (|2.6D . 

Let the impedance matrices seen at the ports of the networks labeled N and Ni in Fig. (1121) be Z and Zi 
respectively. Assuming that Z is positive-real the aim of the resistance and Belevitch transformer extraction 
is to get Zi PR together with Re [Zi.n (jwi)] = 0 for some frequency uji where Zi.n is the (1,1) entry of the 
impedance matrix Zi. As we will see below those are necessary and sufficient conditions for the multiport 
Brune algorithm. 

3.1.2. Resistance Extraction 

The most direct way to get Re [Zi.n (jwi)] = 0 is to first make an eigenvalue decomposition for the 
Hermitian part Z// of Z at each frequency w 
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Zh ijio) = U{uj)Siio)U^ (w) 


( 3 . 10 ) 


where 


Zh ijuj) = -{Zijio) + Z^ i-juj)) 


(3.11) 


S (w) can be assumed to be the diagonal matrix having the eigenvalues aj (a;)’s for 1 < j < iV on its diagonal 
in increasing order (uj (lo) > 0 since Zh (jw) is positive semi-definite): 


/ ai {oj) 


S (w) = 


V 0 


UN (w) 

and U (w) can be chosen to be orthogonal. 

If we then choose the following value for the extracted resistor ri 

ri = min ai (w) 

0<a;<oo 

with wi being the frequency at which the minimum in Eq. (13.131) occurs 

ai (wi) = ri 


We have Re [Zi^n (jwi)] = 0 with 




Zi = (oji) ZU (uJi) - ri 


One then needs to choose the turns-ratio matrix Ti as 


V 0 


0 \ 
0 / 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 


Ti = C/(wi) 

One can also imagine extracting first the resistor and then the Belevitch transformer. The value of the 
resistor ri in that case can be computed using the following formula 


ri = min A(u;)/Aii(w) 

0<a;<oo 

where A (cu) is the determinant and An (w) is the (1,1) minor of Zh (jw)- 

Let Wo be the frequency at which the minimum in Eq. (13.171) occurs such that 


(3.17) 


ri — A (wo)/All (wq) ( 3 . 18 ) 

Then the Belevitch transformer matrix Ti is given by the matrix that simultaneously diagonalizes 
/I 0 \ 


Zff (juJo) and 


0 


V 0 


such that 


0 y 


Zh (/wo) = TiDTf (3.19) 

where D is a diagonal matrix with D (1,1) = ri. Ti in that case can be found using the Gauss diagonalization 
procedure (3^ . 

The formula given in Eq. (13.171) is nice in the sense that it doesn’t require an eigenvalue decomposition 
for each frequency w, 0 < w < oo. However the Belevitch transformer matrix Ti obtained by the above 
method is in general non-orthogonal. 
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3.1.3. Extraction of the Reactive Part of the Multiport Prune Stage 

Let the subnetwork N 2 in Fig. m be described by the following state-space equations 


X2 — A 2 X 2 -f B 2 U 2 (3.20) 

y2 = C 2 X 2 + D 2 U 2 (3-21) 


where 


and 


B 2 = B 


2A 


B 


2B 


C 2 = 


C 2 A 

C 2 B 


D 2 = 


D 2 AA 

D 2 BA 


D 2 AB 

D2_BB 


U2 = 


U2A 

U2B 


y2 = 


y2A 

y2B 


(3.22) 

(3.23) 

(3.24) 

(3.25) 

(3.26) 


where U 2 A is the current into the first port of the subnetwork N 2 in Fig. (fl^ and U 2 _b is the vector holding 
the currents at the remaining ports (ports 2 — N) of the subnetwork N 2 . Similarly y 2 A is the voltage across 
the first port of the subnetwork N 2 and y 2 B is the vector holding the voltages across the remaining ports 
(ports 2 — N) of the subnetwork N 2 . {A 2 , B 2 , C 2 , D 2 } is a state-space realization for the impedance Z 2 (s) 
seen at the ports of the network N 2 . 

Then the network Ni is described by the the following equations 


i 

V iii 


f A 2 
0 

C 2 A 

\ ni \/L\ 


0 

0 

1 

niA/XiCT 


B 2 A 

ni ^/L\C\ 

_ D 2 AA 

n'jLi 


X2 

XCi 

XL^ 


( B2A/ni B 2 B 

1 —l/ni vf 

V Cl V Cl 

O'?. A A D 94 Fi 

\ n\PTE niy/Ll 


UlA \ 
UlB J 


(3.27) 


f Via 

V yiB 


C2A 

l-l/ni 

D 2 AA \ 


ni 

VCT 

ntVil 1 

XCi 

V 

C 2 B 

y/Ci 

ni \/Li / 


from which we identify 


D2AA/n\ 'D2ABlni 
^2BA/ni D 2 BB 


Ai = 


A 2 

0 

C 2 A 


V "^1 


Bi 


Cl = 


/ B2A/'ni 
1 - 1/711 

y/Ci 
. D 2 AA 
\ n\y/^ 


ni v^LiCi 

B 2 B 


B 2 A 

ni y/LiCi 

_ D 2 AA 

niLi 


y/ Cl 
D 2 A R 


C 2 A 

ni 


1/y^i D 2 AA 


y/Cl 

L/1 


Do R A 


UlA 

UlB 


(3.28) 


(3.29) 


(3.30) 


(3.31) 
































and 


/ D2AA/n\ T)2AB/ni \ 
\ ^2BA/ni D 2 BB ) 


(3.32) 


xi = ^ xci j (3.33) 


where xci = \fC\vcYi ^Li = y/LaiLi and Ui = {vi 2 , ■ ■ ■, vin)^uia is the current into the first port of the 
subnetwork Ni in Fig. m and UiB is the vector holding the currents at the remaining ports (ports 2 — N) 
of the subnetwork A^i. Similarly yiA is the voltage across the first port of the subnetwork and yis is the 
vector holding the voltages across the remaining ports (ports 2 — N) of the subnetwork Ni. {Ai, Bi, Ci, Di} 
is then a realization for the impedance Zi (s) seen at the ports of the network Ni. 

Now we will state the multiport version of the fundamental lemma stated in Section (1^ to show how 
to transform state-space equations given for a minimal realization of the impedance Zi into the form in Eqs. 
(|3.27|) and (|3.28|) . For details refer to [l^ . 

The Multiport Synthesis Lemma 

Let {Aa,Ba,Ca,Do} be a minimal realization corresponding to the positive-real impedance Zi (s) satis¬ 
fying Zi^ii (jwo) -I- Zi ii {—juJo) = 0 for some frequency ujo, Zin is the (1,1) entry of the impedance matrix 
Zi (We also assume that jwq is not an eigenvalue of A^). Then there exists a coordinate transformation 
matrix T such that Ai = TA^T'^, Bi = TB^, Ci = and Di = are of the form given in Eqs. 


To compute T we will follow the algorithm described in Fundamental Lemma in Section (12.51) . Before 
applying the algorithm we set 


ba = B^ei (3.34) 

cl = Cle, (3.35) 

and apply the one-port algorithm described in the Fundamental Lemma in Section (12.51) to the set {A^, ba, Cq, Da} 
where ei = ( 1 0 ... 0 ) ; that is we apply the one-port algorithm by picking up the first columns of 

Ba and Cl matrices. 

To see why {AiJBi, Ci,Di} is an equivalent realization for the impedance Zi (s) see Section (12.21) or 
Theorem (3.3.9) in [l^ . 

3.1.4- The Multiport Capacitive Degenerate Stage 

Similar to our discussion in Section (12.5.11) it is possible that the minimum in Eq. (13.141) occurs at 
infinity, uji = 00 . Such a case needs the extraction of the degenerate reactive stage shown in Fig. (1151) . As 
in the one-port case the reactive Brune stage doesn’t involve any inductive part. Note also that we don’t 
have the iz-type transformer coupling the first port to the remaining ports. 

To synthesize such a stage we need to modify our treatment in Section (13.1.31) . Assuming that the 
subnetwork N 2 in Fig. is described by Eqs. p.20l) - (l3.26l) . the network Ni is described by 


*2 

ici 


f Via 

V yiB 


A 2 — 


1 C 2 A 


D 2 AA 
C 2 A 


D 2 AAVC 1 


0 

/-I ID^raC^a 

- dTIT- 


B 2 A \ 

D 2 AAy/Ci j 

( 

D2AAC1 J 

V , 

1 \ 



1 ( 

D2RA J 

^2AaVCT / 

^ V 
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B 2 R — 


^2Ai->2AB 


D 2 AA 

Dl A B 

D2AAy/Ci 


UiA \ 
UiB J 



0 

r>2BB - 


UlA \ 
UiB J 


(3.36) 


(3.37) 



























Figure 15: Extraction of the multiport capacitive denegenerate stage in the multiport Brune circuit when wi = c» in Eq. 
We observe that the i/-type transformer is not necessary in that case. 


from which we identify 


Ai = 


B 2 AC 2 A 


C-2A 


D 2 AA 


D2AAVC1 


B 2 A 

D2AA'/Ci 

1 

1^2AaCi 


Bi = 



B 


2B 


_ B2AD2AB 
D 2 AA 
Dl A B 


D2AA s/Cx 


Cl = 


^ _ D2BAC2A D2BA 

D2AA 


D 2 AAVC 1 


Di = 


0 

0 D 


2BB 


0 

D2BaD2>^ B 

D 2 AA 


and 


(3.38) 

(3.39) 

(3.40) 

(3.41) 


Xl = 


X2 

XCi 


(3.42) 


where xci = \/Civci\ uia is the current into the first port of the subnetwork Ni in Fig. m and UiB is 
the vector holding the currents at the remaining ports (ports 2 — N) of the subnetwork Ni. Similarly yiA 
is the voltage across the first port of the subnetwork iVi and yiB is the vector holding the voltages across 
the remaining ports (ports 2 — N) of the subnetwork Ni. {Ai, Bi, Ci, Di} is then a realization for the 
impedance Zi (s) seen at the ports of the network iVi. 

One needs to modify also the Multiport Synthesis Lemma as follows: 

The Multiport Synthesis Lemma (multiport capacitive degenerate case) 
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Figure 16: Extraction of the multiport inductive degenerate stage. 


Let {Aa,Ba,Ca,Da} be a minimal realization corresponding to the positive-real impedance Zi (s) sat¬ 
isfying Zi ii (jwo) -I- Zi^ii {—juio) = 0 for some luq = oo, Zi^n is the (1,1) entry of the impedance matrix Zi. 
Then there exists a coordinate transformation matrix T such that Ai = TA^T^^, Bi = TB^, Ci = CaT~^ 
and Di = Dq are of the form given in Eqs. 1 (I3.38I) - (I3.411) 1. 

To compute T we will follow the algorithm described in The Fundamental Lemma (one-port capacitive 
degenerate case) in Section (12.5. 1|) . Before applying the algorithm we set 

ba = B^ei (3.43) 

cl = Clei (3.44) 

and apply the one-port algorithm described in The Fundamental Lemma (one-port capacitive degenerate 
case) in Section (12.5. 1|) to the set {A^, b^, c^, Da} where ei = ( 1 0 ... 0 ) ; that is we apply the 

one-port algorithm by picking up the first columns of Ba and Ca matrices. 

3.1.5. The Multiport Inductive Degenerate Stage 

Similar to our discussion in Section (|2.5.2I1 for the one-port Brune circuit extraction there is the possibility 
of the minimum in Eq. (13.1411 occuring at wi = 0. Such a case corresponds to the extraction of an inductive 
degenerate stage as shown in Eig. dm). Again it is straightforward to extend the multiport state-space 
Brune algorithm presented in Section (13.1.31) to extract the multiport inductive degenerate Brune stage in 
Fig. (dH); we do not describe this algorithm here. We would like to also note here that one can extend the 
treatment in Section (12.91) of the one-port inductive degenerate case for the quantization and dissipation 
analysis to the multiport case in a straightforward way provided that the resistance ri in Fig. © is high 
(that is the loss introduced by the resistor ri is small). 

3.2. Quantization of the Multiport Brune Circuit 

The multiport Brune circuit contains ideal transformers. In this section we show that one can eliminate 
transformer branch variables and write a set of effective Kirchhoff relations for the rest of the branch currents 
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Figure 17: Multiport Brune stage. Tree branches are shown in black and chord branches are shown in blue. Note that the 
Belevitch transformers Tj and Tj+i are reflected compared to the Fig. JTSj. 


and voltages. The effective Kirchhoff relations are given by the loop matrix involving turn ratios which 
we define below. The treatment here is similar to the analysis done in Section (12.61) for the one-port state- 
space Brune circuit. However as we will see the addition of Belevitch transformers and u-type transformers 
makes the analysis more involved. We note that we replace the shunt resistors Rj’s in the last stage by 
capacitors Cr^ ’s for 1 < j < iV as shown in Fig. ((^(11) to do the dissipation analysis as discussed in detail in 
Appendix (I7.1.2L 

We will follow the same approach of Section (|2.6I1 . We write again the Kirchhoff’s laws (Eqs. (17.51) . 
(17.61) 1 for the multiport Brune circuit in Fig. m whose stage is shown in detail in Fig. (da 

FIch = -Itr (3.45) 

F^Vt, = Veh (3.46) 

where we have again assumed that there are no external fluxes in circuit loops. As in the one-port case F 
is again the loop matrix with entries being 0, 1 or —1 derived by a graph theoretical analysis of the circuit 
Q. After the effective Kirchhoff analysis done below F will be replaced by the effective loop matrix 
with real-valued entries. 

Itr and Ich are the tree and chord branch current vectors in Fig. (fTTl) respectively partitioned as follows 


Itr= (lj,lL,Iz,lf^) (3.47) 

Ieh= (3.48) 

and tree and chord branches’ voltages are partitioned respectively as 

Vtr= (3.49) 

Vch= (3.50) 
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Here labels J, L, Z^ C, T correspond to Josephson junction, inductor, resistor, capacitor and ideal trans¬ 
former branches, respectively. 

Our aim here is to write an effective set of Kirchhoff relations as 

fO/iO/ = _iO/ (3.51) 

(pT) O/ ^e// ^ ^e// (3 g2) 

where transformer branches are eliminated such that 


■e// 

tr 

= (Ij,Il,Iz) 

(3.53) 

■eff 

'ch 

= Ic 

(3.54) 

■eff 

tr 

= (Vj,Vi,Vz) 

(3.55) 

■eff 

ch 

= Vc 

(3.56) 


We note here that the entries of the effective loop matrix in Eq. (13.511) are real numbers (as in the 
one-port case) being functions of ideal transformer turn ratios as we will see below. 

Here we will do this effective loop matrix analysis for the Kirchhoff’s current law to get the matrix 
in Eq. (13.511) . It is important to note that one should also do a similar analysis for the Kirchhoff’s voltage 
law to get an effective in Eq. (13.521) and verify that 


(fT)^// ^ (F®//)'^ (3.57) 

holds. This we do in the Appendix (17.2.21) . Eq. (13.571) is important to keep the various matrices of interest 
like the capacitance C and stiffness Mq matrices symmetric. 

To show how one can find such a F®-^-^ matrix we will further partition transformer current vectors 
and in Eqs. (13.471) and ()3.48|) . We first note that left branches of all transformers in the circuit in Fig. 
(Ha are chord branches (colored in blue) and that right branches of all transformers are in the tree (shown 
in black in Fig. mi)). Hence we can write 


where 


and 


t(*^) _ 

Lrp — 

(l.'/'t 

4'”, 4®’) 

(3.58) 

II 

1—1 


4>,i<,«) 

(3.59) 

t(-R) — 

(«',. 

...47) 

(3.60) 

t(-^) _ 

Irj, — 

(47,. 

■■.47„) 

(3.61) 

t(-R) — 

u 

(47., 

...47) 

(3.62) 

T(^) - 

('17. ■ 

...47) 

(3.63) 

riL) _ 
— 

(47,. 

...47«) 

(3.64) 

t(^) - 

(47.. 

...47) 

(3.65) 
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Figure 18: i^-type transformer circuit at multiport Brune stage in Fig. i(T7) i. We also show the connection of the transformer 
to the capacitor Cj to make the correpondence to the Fig. (Eli more clear. 


with 


j(L)iR) 


t(«) 

1^7 


/ 

V 

I 

V 


i: 


{L){R) 


T.-.l 


i: 


iL){R) 


Tj,N 
l(R) \ 


jiR) 


(3.66) 

(3.67) 


where are vectors of length for 1 < j < M + \ and are vectors of length (iV — 1) for 1 < j < M. 

Before moving further in the analysis we beriefly review the relations between the currents through left 
and right branches of the three different types of transformers in the multiport Brune circuit. We also give 
voltage relations for completeness although we don’t need them for the analysis in this section. However we 
will refer to the voltage relations in the Appendix (17.2.211 . 

ji^)(R) .g vector of currents through the left(right) branches of the Belevitch multiport transformer 
appearing at the multiport Brune stage, 1 < j < M + 1 as shown in Fig. dni). Hence by Eqs. (13.81) and 
(1^ we have 


(3.68) 

(3.69) 


where is the vector of voltages across the left (right) branches of the Belevitch multiport transformer 

appearing at the multiport Brune stage, 1 < j < M + 1 in Fig. IHD- 

The turns ratio vector of the iT'-type transformer i/j at the multiport Brune stage is given by 


^i.2 


V ^j,N 


(3.70) 


The detailed circuit representation of the i^-type transformer at the stage is given in Fig. (fTSl) . The 
current on the left branch is related to the currents on the right branches by 
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Figure 19: Conventions for the direction of currents for the n-type transformer of stage in the multiport Brune circuit in 
Fig. {TtJ- The inductive branch Lj is shown for orientation purposes with Fig. 


/ 




(3.71) 


and the voltages across the right branches of the i^-type transformer are related to the voltage across its left 
branch by the following formula 


( \ 


(R) 




(3.72) 


\ X.Af / 


The detailed circuit diagram with current direction conventions for the n-type transformer is shown in 
Fig. (dH). For this type of transformer the relations between currents and voltages are given by 


(3.73) 

(3.74) 

Now we can begin our analysis. We will proceed from the last (rightmost) stage to the first(leftmost) in 
Fig. (dl]) by relating the currents of consecutive stages. 

Starting at the last stage we have (see Fig. dnn) 

Itm+1 = (3.75) 

where Icr are the currents through the capacitors (substituted for the shunt resistors) of the last stage. 

The currents of inter-stage transformers are given by Eq. (13.681) 


4";) = (3.76) 

for 1 < j < M -f 1, where T^- is the {N x N) Belevitch transformer matrix of the stage. The currents of 
consecutive inter-stage transformers are related by 


for 1 < j < M, where 


t(^) 


-ei/c, 


J ij + i 


/ 1 



• • ■ 

1 0 
0 


1 / 


(3.77) 


(3.78) 
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Figure 20: Last stage of the multiport Brune circuit. Tree branches are shown in black and chord branches are shown in blue. 
Shunt resistors Rj's in the last stage are replaced by capacitors C/^.’s for 1 < J < • Their contribution to the dissipation is 

computed in Appendix 117. 1.21 1. 


Aj is a {N X N) matrix and ei is the unit vector ei = (1, 0,..., 0)^ of length N. We note that in the case 
of a degenerate stage as in Fig. m for the stage we have = 0 for 2 < fc < iV , hence Aj is the 

identity matrix in such a case. 

The current II^ through the inductor Lj at the multiport Brune stage can be written 

lL,=-Ic,+v,lir% (3.79) 

where Vj is a row vector of length N : 

Vj = {{l-nj) -Vj2 ■■■ -Vj,N ) (3.80) 

Using the Eqs. (13.751) . (|3.76l) . (|3.77l) and (13.791) we can iterate over the index j backwards starting at 
j = M + 1: 


I 


{L) 

Tm+1 

^Lm 


1 M 
1 M 


= -Tm+iIcji 

-Icm + = -^Cm - VmTm+iICr 

-ei/cM + AmIt2+i = - AmTm+iIck 

TmIt^ = —TmGi/cm ~ TmAmTm+iIch 
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^Lm-1 
t(^) 


Tm-1 

(L) 

Tm-1 




= —Icm-1 + “ ~Icm-i ~ I’M-iT^MeilcM ~ '<^m-iTmAmTm+iIch 

= —Gi/Cm-i + “ ~®1^Cm-i ~ Am-iTmGi/cm — Am-iTmAmTm+iIcji 

_ rp t(-^) 

= —Tm-iGi/cm-i ~ Tm-iAm-iTmGi/cm — Tm-iAm-iTmAmTm+iIck 
= -^Cm-2 + 

= —ICm-2 ~ 'I^M-2Tm-1Gi/c'm-i ~ '*^M-2Tm-iAm-iTmGi/cm + 

—i^m-2Tm-iAm-iTmAmTm+iIcr 


(3.81) 


Hence we conclude that one can write 


II = -F^"lc 


(3.82) 


with for 1 < j < M : 


pie// 

pie// 

pie// 

pie// 

pie// 

pie// 


(j, k) = 0 

U, k) = 1 

(j,fc) = T-jTj+iei 

(/) fc) = ... Tfe-iAfc-iTfcGi 

{j, k) = UjTj+i Aj+i... TmAmTm+iG^-m 
(/, k) = vmT^m+i^u-m 


for k < j 
for k = j 

for k = j + 1, and j < M 

for j + 2 < k < M 

for j < M and M + l<k<M + N 

for j = M and M + l<k<M + N 


(3.83) 


where efc is the unit vector of length N non-zero only at its entry such that {j) = 0 for j k and 
Gfc (fc) = 1. We assumed the following ordering for the capacitors 

{Cl,..., Cm, Ciji,..., (3.84) 

To compute note the following 


Ir. = L 


(fl) 


W.i 

- 

— e;^ 


for 1 < j < M. 

Referring back to the iteration in Eqs. (13.811) and using Eq. (13.861) we can write 


Iz = -F|"lc 


with for 1 < j < M : 


(F%ff{j,k)=0 
lF^"(j,fc) = l 

I ^‘zc Oi k) = e{’AjTj_|_i ... Tfe_iAfc_iTfcei 

Oi k) = ef AjTj_|_i ... TMAMTM-i-iGfc_M 


for k < j 
for k = j 

for j + 1 < k < M 
for M + l<k<M + N 


(3.85) 

(3.86) 


(3.87) 


(3.88) 
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Figure 21: To compute for simplicity, we assume that all the ports of the multiport Brune circuit are shunted by 

Josephson junctions Lj ^, • • •, • 


To compute we will assume for simplicity that all the ports are shunted by Josephson junctions as 
shown in Fig. m so that all the junctions belong to the spanning tree. Later we will allow for resistors 
and voltage sources shunting the ports of the multiport Brune network. We note the following 

1/ = (3.89) 

Hence referring back to the iteration in Eq. (I3.81|) we deduce using Eq. (I3.89|) 


Ij 


tH) 

^Ti 


—Tiei/ci 


TiAiT2ei/c2 ... — TiAi .. .Tm-iA.m-i'Tm^iIcm + 


—TiAi... TmAmTm+iIch 

Hence for 1 < j < A 

for k = 1 
for 1 < k < M 
for M + l<k<M + N 


'F%f {j,k) = ejT,e, 

‘ F"/c ih k) = ejTi Ai... Tfc_i Afc_iTfeei 
Tjc/ (.?’ k) = ejTiAi... TMAMTjvf+iefc-M 


where Fj-^ is defined by 


(3.90) 


(3.91) 


I; = -F}"lc (3.92) 

Now that we derived effective loop matrices we will follow the Appendix (17.1.11) and (17.1.21) to compute 
Hamiltonian matrices and to do dissipation analysis due to the resistors in the multiport Brune circuit. We 
will repeat here some of the definitions - which will be used with the effective loop matrices- of the Appendix 
(17.1.11) for convenience. 

With the effective loop matrix 


Tre// _ 

■^C — 


■pe// 

JC 

pe// 

LC 


(3.93) 
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we can compute the capacitance matrix C defined in Eq. (17.131) 


C = 


Cj 

0 


0 

0 





(3.94) 


where we have the following partitioning to identify the submatrices Co and Cr corresponding to ordinary 
capacitances and capacitances to be replaced by shunt resistors in the last stage of the multiport Brune 
circuit, respectively 


with 


0° C* 


Cl 


and 


Cr = 


V 0 

/ Cr, 

V 0 


Cm 


0 \ 


Cr^ 


and the matrix Cj in Eq. (13.941) for the Josephson junction capacitances is 

/ Cj, 0 


C,7 = 


V 0 


Cj. 


Partitioning also according to the partitioning in Eq. (13.951) as in Eq. (17.391) 


= ( JcCf Tceff 


Co -^Cr 


we can decompose C in Eq. (|3.94|) as in Eq. (17.421) 


C = 


Cj 0 
0 0 


vs"c 




— Co + Cr 

where we defined as we did in Eqs. (17.431) and (17.441) 


(3.95) 

(3.96) 


(3.97) 


(3.98) 

(3.99) 

(3.100) 

(3.101) 


(3.102) 

Cr = (3.103) 


Co is the capacitance matrix that appears in the system Hamiltonian in Eq. (13.1081) below for the multiport 
Brune circuit whereas Cr is a dissipative term which we will treat later below by making the substitution 
as described in Eq. (17.3711 for the shunt resistors in the last stage, 
and Q matrices defined in Eqs. (|7.33|) . (I7.34|l are simply 
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/ 1/Ll 0 \ 

V 0 1/Lm / 



since there are no chord inductors in the circuit. Hence we have by Eq. (|7.31D 


(3.104) 

(3.105) 


Mo = 


where 0,/ is a Nj x Nj zero matrix. Here since we assumed that all ports are shunted by Josephson junctions 
Nj = N. 

Using the Eq. (17.941) we have the following Hamiltonian for the multiport Brune circuit 

ns = + u ($) (3.108) 

where 




( 


V 0 


1/Li 


0 \ 


I/Um / 


(3.106) 

(3.107) 


[/(^) 



cos {ifj) + 


To treat the dissipation we first note that the diagonal tree impedance matrix 


(3.109) 


Z = 


( ri 

V 0 



(3.110) 


consists of in-series resistances ri,..., tm of the multiport Brune circuit in Fig. (ED. However, as we noted 
in Appendix (17.1.211 we will treat each in-series resistor rj separately for 1 < j < M. So instead of using the 
full matrix dehned in Eq. (13.881) we will use its rows corresponding to the in-series resistors 

Tj’s for 1 < / < M, in our dissipation treatment below. 

We compute the coupling thj of the bath due to the resistor rj for 1 < / < M to the circuit degrees of 
freedom using Eq. (17.681) with the effective row vector F^J^^ (since we treat resistors one at a time F^^/q^ 
is the row of the matrix F^^^ defined in Eq. (17.591) corresponding to the resistor rj) 


and using Eq. (17.801) 


■^Co '^0 


myj 

myL 


(U"J 


CJ = —lUJVi 


1 -h iu:rjFl^/co^o 


{<k) 


which is related to the kernel as derived in Eq. (17.781) (we note that Cz,rj (w) is a scalar) 


(3.111) 

(3.112) 

(3.113) 
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(3.114) 


Kj (w) = -Oj'^Cz,rj (w) 

and the spectrum of the bath due to the resistor rj is given by Eq. (17.811) 


Jj (w) 


Im [Kj (w)] 


UJ^Tj 


1 + 


f C' 


«ky 


(3.115) 

(3.116) 


The contribution of the resistor rj to the relaxation rate is given by the formula in Eq. (17.821) 

^ 1(0 |mj • ^1 1)1^ Jj (woi) coth (3.117) 

where wqi is the qubit frequency. In this multiqubit case, this formula can be used to get to the relaxation 
rate of any qubit, by choosing the appropriate initial and final quantum states. For example, the relaxation 
rate of qubit 1 may be calculated by taking |0) to be the global ground state in which each qubit is 0, and 
|1) to be the state in which qubit 1 is in its upper eigenstate, while the other qubits remain in their ground 
state; we would write this more fully as the state |1000...). Taking this to be the state |0100...) will give Ti 
for the second qubit, and so forth. In fact, this formula can be used to get the relaxation rate between any 
two multi-qubit eigenstates. A further refinement of these prescriptions may be needed if the eigenstates of 
the multiqubit system are entangled. 

To compute the dissipative contribution of the resistors {Ri,..., R]\j} shunting the last stage we make 
the substitution C/j ^ (*w)~ done in Eq. (17.371) to get the resistance matrix defined in Eq. (|7.48l) 


where R is defined as 




Tre// 


R 


-1 



/ Ri 0 \ 

\ 0 Rn / 


(3.118) 


(3.119) 


We again treat the shunt resistors {i?i,..., Rn} one at a time, that is, to compute the contribution of 
the resistor Rj to the relaxation rate we set = oo for 1 < and k j and we short circuit in series 

resistors by setting rfc = 0 for 1 < fc < M. The coupling vector which couples the bath due to Rj to 
the circuit degrees of freedom is given in Eq. (17.541) as 


mfi, = r: 


■eff 

Rj,CR 


(3.120) 


where R^I.^Cr column of the matrix and the dissipation kernel due to Rj is given in Eq. 

(I73HD as " "" 


KRj ( w ) 


ioj 

Rj 


and the spectral density Jr^ of the bath due to Rj is given by Eq. (17.561) 


(3.121) 


JRj {^) = Im [Kr. (w)] 


ut 
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(3.122) 

(3.123) 
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Figure 22: The multiport Brune circuit shunted by a resistor at its port and by Josephson junctions at the remaining ports. 

Hence by Eq. (|7.58p we compute the contribution of the resistor Rj to the relaxation rate as 

= 11{0 |m«, - «| 1)1’ Ju, Mcoth (3.124) 

3.3. Resistors Shunting the Ports 

One can also think of shunting some of the ports in Fig. m by resistors. For simplicity let’s assume 
that we shunt only the port of the multiport Brune circuit by a resistor and the rest of the ports 
by Josephson junctions as shown in Fig. (EH). The branch of the resistor Rj belongs to the spanning tree 
hence its treatment will be similar to the treatment of the in-series resistors rj in the previous section. We 
will derive an effective loop row-vector F^^^ corresponding to Rj. Before doing this we note that the 

row of the effective loop matrix F^/J derived in Eq. (13.911) above needs to be dropped since we replaced the 
X Josephson junction in Fig. EID by the resistor i?J. 

We first note the following 


r _ AL) 
- -'Tij 


Hence from the iteration in Eq. (I3.81|) we conclude 


(3.125) 


% 


—ejTiei/ci — ejTiAiT2ei/c’2 ... — ejTiAi ... Tm-i AM-iTMei/cM + 
— ejTi Ai... Tm AmTm-i-iIcr 


From which we conclude that 


(3.126) 


— F^/pIq 


where F^^J is a row vector of length (M -|- N) with 
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(3.127) 


(^) = eJ’Tiei for k = 1 

* ^r/,C (^) = ejTiAi... Tfe_iAfe_iTfeei for 1 < k < M 

^ G (^) “ gTTi Ai... T M ■^m'F M+i^k-M f or M + l<k<M + N 

Vi’ 

So one then needs to append F^j^x^q to F^Jq defined in Eq. (13.881) as its last row. Also Rj will appear at 
the last diagonal entry of the tree impedance matrix Z defined in Eq. (13.11011 


t ri 

0 \ 

V 0 

tm 

R- y 


(3.128) 


To compute the contribution of i?| to the dissipation rate we follow the same procedure done for the 
resistors r^, 1 < j < M in Eqs. (I3.111l) - (|3.117ll above and compute the coupling of the bath due to 

the resistor to the circuit degrees of freedom using Eq. (17.6811 with the effective matrix F^Jx^q^ 




_ TT^//O 

- ^Co '-^0 




^RfJ 

mR-L 


where we assume the following partitioning for 


pie// _ / pie// 

^ R^,C — ^ 


fie// 


where F^j^x^^^ and are row vectors of length M and N, respectively. 

Using Eq. (17.801) 


(3.129) 

(3.130) 


(3.131) 


Cz.R^ (^) = 


1 + iuR^jF 


e// 

RfCo 


Co 



-1 


which is related to the kernel as derived in Eq. (17.7811 


Krx (w) = —uj'^Cz,R^ (w) 

and the spectrum of the bath due to the resistor Rj is given by Eq. (17.811) 


(3.132) 


(3.133) 


Jflx (w) 


Im 


Kr^ (w) 




l+a;2 {R^y 


f 

flTCo'^O 



And by Eq. (17.8211 the contribution of R^ to the loss rate 


(3.134) 

(3.135) 


1 




4 

h 


mflx • # 


2 

Jr^ (wol) coth 


hujol \ 
2kBT) 


If we have several resistors shunting the multiport Brune circuit the analysis is straightforward. One then 
needs to drop the rows in the F^J matrix corresponding to the ports shunted by resistors and repeat the 
dissipation analysis above for each R^. A full F^J/ is obtained by appending the rows F^x'^^ corresponding 
to each which are appended to Z. 
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Figure 23: Voltage source Vs (t) in series with a resistance Rs shunting the port of the multiport Brune circuit. 


3.4- Voltage Sources Shunting the Ports 

Now we consider shunting one of the ports by a possibly time-dependent voltage source 14 (t) as shown 
in Fig. ((^51) . We also assume that the voltage source has an in series resistance Rg. Assuming also that 14 
and Rs shunt the port the treatment of the resistor Rg follows the same analysis we did above for the 
resistor in the previous section. To treat the voltage source we first note the following 

/v, =/g^. (3.136) 

Hence with the same loop analysis we did for the resistor above we can derive an effective loop matrix 
Fy/c such that 

ly^ = (3.137) 

where with F^-^^ being defined in Eq. (|3.127L 

Now Fy^^Q will appear as one of the rows of the Fyc matrix. One can then use Fyc to first compute 
my defined in Eq. (17.201) and then the coupling vectors Cy and Cy defined in Eqs. (17.141) and (17.251) . 
respectively, in Appendix (17.1.111 . Making the substitution C/j v- (ia;)~^R“^ in Eq. (17.371) and following 
the analysis at the end of Appendix (17.1.11) we get the following Hamiltonian derived in Eq. (17.941) for the 
multiport Brune circuit shunted by the voltage source at one of its ports 

ns = \{Cl- (Cyo + CyR (t)) * Vy {t)f (Q - (Cyo + CyR (t)) *Vy{t)) + U ($) (3.138) 

where the voltage source vector Vy (t) has 14 {t) at its corresponding entry. The matrices Cyo and Cy^R {t) 
are defined, respectively, in Eqs. (17.851) and (17.9211 in Appendix (|7.1.3I) which we repeat here for convenience 

Cyo=T^J/Co(F‘y%y (3.139) 

Cv,R {^) = Cv,R {^) + Cy (uj) (3.140) 
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Figure 24: 2-port 1-stage multiport Brune circuit. The part of the circuit labeled N contained in the dotted box is the 2-port 
Brune circuit shunted by two Josephson junctions at each of its ports. 


where in the definition of the frequency-independent coupling matrix Cy^o we assumed the following parti- 
tioning for Fy^ as in Eq. (17.831) 

Fvc = ( ) (3-141) 

For the definitions of the matrices Cv,r (uj) and Cv (w) appearing in Eq. (13.1401) we refer the reader to Eqs. 
(17.861) and (17.911) . respectively, in Appendix (|7.1.3I) . We recall that effective loop matrices should be used in 
those definitions. 

We would like to note here that the analysis above is done only to obtain the time-dependent Hamiltonian 
in the case of a AC voltage drive 14 (t ); we will only use our formalism to compute Ti times only for unexcited 
systems, i.e. when there is no time-dependent voltage applied. 


4. Examples 

Here we illustrate the method described in the previous section with a generic 2-port 1-stage Brune 
circuit. 

4-F 2-port 1-stage Generic Brune Circuit 

We consider a generic 2-port 1-stage Brune circuit shown in Fig. (1^ . The 2-port Brune circuit is shown 
in the dotted box. This 2-port Brune circuit is shunted by two Josephson junctions Lj^ and Lj^ at each 
of its ports. There is only one reactive stage which corresponds to a single pole at a finite frequency in the 
partial fraction expansion for the impedance fit in Eq. (HD. 

First we need to replace the resistors Ri and i ?2 in the last stage by capacitors Cr^ and Cr^ as shown 
in Fig. (1^51) . We will later make the substitutions Cr^ ^ ^ii 2 ^ 77^7 dissipation analysis. 

There are three capacitors in the circuit in Fig. (I^Hl) . All the capacitors are in chord branches. The 
currents through those capacitors are given by the vector 

Ic = f iZ, ] (4.1) 
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Figure 25: Modified 2-port Brune circuit. Chord branches are shown in blue, tree branches are shown in black. 


As described in the previous section we will eliminate transformer branches and write the cnrrents through 
tree branches as functions of capacitor currents. This way we will compute effective loop matrices F^/J, 
Flff and . 

As we see in Fig. (1^51) the two Belevitch transformers Ti and T 2 are given by the turns ratios matrices 


We can now write 


riL) 

^2 


Ti = 
T 2 = 



( ^11 

tl 2 

\ ^21 

t 22 

f Ull 

U12 

V U21 

U22 

II 

-T2IC 

Mil 

U12 \ 

U2I 

U22 ) 


Icr^ 


where we used I 


(R) 

T 2 


—Icr- The current through the inductor Ti is given by 


(4.2) 

(4.3) 


(4.4) 

(4.5) 


= 


from which we identify 


-Ic, + 


-/Ci - ((1 - Ui) Uii - U 21 V 12 , (1 - «l) U 12 - U 22 V 12 ) 


Icr, \ 
^Cr^ ) 


(1 “ “11 “ “211^12, (1 - ni) U 12 - U 22 V 12 ) 

where F^-^ is defined by the relation 


(4.6) 

(4.7) 


(4.8) 
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(4.9) 


h, = -F^"lc 


We now move to the transformer Ti and write 


r(«) _ 


-Ic,ei + 


-Ic, 

-ICl 


-Vl2 

1 

-1^12 

1 


r(i) 


Mil Ul2 
U2l U22 


ICr^ 

^Cr, 


and 


t(L) 


rp TCfl) 

lllxi 


= -Ic 


—-^^CiTiei - 

til 

t21 


TiAiI^t^J 


til 

t21 


tl 2 

t22 


= -ICl 


til 

t21 


1 -1^12 

0 1 


Mil 

U21 


Ul2 

U22 


tllMii + M 21 (ti 2 — tiil^i 2 ) tiiMi 2 + U22 (tl 2 — tiil/i 2 ) 
t21?^ll + M2I (t22 — t2lJ^12) t2lMi2 + U 22 (t22 — t2lI^12) 

Using the relation Ij = and the Eqs. (I4.13|I - (I4.15|) above we identify 


Icr^ 

^CR2 


Icr^ 

Icr^ 


TT®// _ 

^ JC — 


til tiiMii + M21 (tl2 — tlll^l2) tiiMi2 + M22 (tl2 — tlll^l2) 
t21 t2lMii + M21 (t22 — t2lJ^12) t2lMi2 + M22 (t22 — t2lJ^12) 


where Fj-^ is defined by the relation Ij = —F^/Jlc- 

Using Eqs. (I4.8|l . (14.161) and making the partitioning in Eq. (13.991) we get 


ire// _ 

•^Co - 



and hence by Eq. (13.1021) 


Co = 


Cj 

0 


■^Co '-^0 




Cji+tfiCi tiit2iC'i tiiCi 

t2ltllC'i Ujj+tliCi t2lCi 


tiiCi 


t2lC'i 


Cl 


where we noted that Cq = Ci and assumed that Cj = 

necessary here to have a non-singular Cq matrix. 

Mo is given by Eq. (I3.107|) as 


Cji 

0 


0 

Cj, 


. We also note that a non- 


Mn = 



0 

0 

1 

Li 


(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 

zero Cj is 

(4.20) 














Hence the Hamiltonian is given by Eq. (|3.108p as 


ns = ^Q'^Co^Q + U{^) (4.21) 

with 

C/W = -(^) L7icos(vP^) + i$^Mo$ (4.22) 

and 




JNotinsf also R. = we identifv bv Eq. (14.121) 


1 “11 “ “211^12 M12 - U22’^12 ) 

(4.24) 

is defined by 1^ = —F^^Ic- We note 

K^c = 1 

ri,Co 

(4.25) 

and usine Eq. p.lll|) and Eqs. (|4.17p. (14.251) we comoute 


mi = 


(4.26) 

= 

(‘r) 

(4.27) 

Usiner Eq. (|3.113|) and (|4.25|) we comoute 



Cz,ri (w) = -iojri 


-1 

(4.28) 

= —iujri [1 + iujriCi] ^ 

(4.29) 


and by Eq. (13.11411 we have 


A'i(w) = -io^Cz,rA^) (4.30) 

= iAri [1 + iujriCi] ^ (4-31) 

hence by Eq. (13.1161) we get the spectral density of the bath due to the resistor ri as 


Ji (w) = Im [Ki (w)] 

riA 


(4.32) 

(4.33) 


The contribution of ri to the loss rate is computed then by using the formula in Eq. (13.11711 . 

To do the dissipation analysis for the last two shunt resistors Ri and i ?2 we first note the following, using 
Eqs. (14.811 . (14.161) and making the partitioning in Eq. (|3.99|1 
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tiiUii + U21 {tl2 — tllVl2) tllUl2 + U22 {tl2 — illl^l2) 
~ I ^21^11+^21(^22—^21^^12) t2lUl2 + U22 {^22 — ^21^12) 


(1 - ni)Mii - U2lJ^12 


(1 - ni)Mi2 - U 22 V 12 


(4.34) 


■eff 


In Eq. (I3.12()|) is defined as the first column of and is defined as the second column of T^Cr 
so that 


and 


illMll + U 2 I {tl2 — tllVl2) 
m_Rj = I t2lMll + M 2 I {t22 - t2lVl2) 
(1 - ni)uii - M21I^12 


fllMl 2 + U22 {tl2 — tllh'12) 
^R2 = 1 ^21M12 + U 22 (^22 — t2lVl2) 

(1 - ni)ui2 - U22V12 


(4.35) 


(4.36) 


The spectral densities Jr^ and Jr^ of the baths due to the resistors Ri and i? 2 , respectively, are given 
by the formula in Eq. (13.1231) as 


Jri {^) = 


UJ 


(4.37) 


and 


JR2 (^) = 


UJ 

R 2 


(4.38) 


Contributions of the shunt resistors Ri and R 2 to the loss rate are then computed using the formula in 

Eq. inmi) . 


4 . 2 . 3-port Data for the 3D-Transmon 

In this section we do a multiport study of the numerical example that has been studied with the one-port 
classical Brune analysis in [ij. Below we repeat the description of the system for completeness. 

To show the application of the synthesis method we have just described, we analyse a dataset produced to 
analyse a 3D transmon similar to the one reported in an experiment at IBM @ . Our modeling is performed 
using the finite-element electromagnetics simulator HESS [l^- Since the systems we want to model admit 
very small loss Blli, they are very close to the border which separates passive systems from active ones. 
Therefore it is necessary to take care that the simulation resolution is high enough to ensure the passivity 
of the simulated impedance. Otherwise the fitted impedance Z (s) does not satisfy the PR conditions in 
Section (EH) meaning that there is no passive physical network corresponding to Z (s). 

The simulated device is a 3D transmon, inserted with appropriate antenna structures into the middle 
of a rectangular superconducting (aluminium) box cavity, which is standard in several labs presently for 
high-coherence qubit experiments. Fig. (1^^ shows a perspective rendering of the device, and Fig. (1271) 
shows an intensity map of the fundamental mode of the cavity. The simulation includes two coaxial ports 
entering the body of the cavity symmetrically on either side of the qubit. HESS is used to calculate the 
device’s three-port S matrix over a wide frequency range, from 3.0 to 15.0 GHz. The three ports are those 
defined by the two coaxial connectors and the qubit terminal pair. That is, the metal defining the Josephson 
junction itself is absent from the simulation, so that its very small capacitance and (nonlinear) inductance 
can be added back later as a discrete element as in Fig. ®. The conversion from the S matrix to Tigim is 
calculated using standard formulas . Here we don’t shunt the coaxial ports by 50D terminations as 

we did in [l|. Rather, we keep the simulated impedance matrix Z^im with 3 ports. 

To obtain the impedance matrix Z(s) as in Eq. (jl.ll) fitted to the numerical impedance matrix Zsim, 
we use the MATLAB package Vector Fitting (VF) [l^. Applying the compacting technique described in 
we get a state-space description for Z (s) of dimension 21 (the dimension before compacting being 51). 
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Figure 26: Geometry of the 3D transmon qubit simulated in HFSS. Light blue is perfect conductor and dark blue is the 
vacuum. The qubit port terminals are defined on a dielectric substrate located at the position of the red line. Two coaxial 
ports are positioned symmetrically on each side of the substrate. The cavity dimensions are (height, length, width) = 
(4.2mm, 24.5mm, 42mm). 


Applying the multiport Bruiie algorithm described in Section (j3.ip we obtain the circuit parameter values 
listed in Table © with the number of stages M = 12. There are three degenerate stages, namely stages 
1,11,12. For the multiport Brune circuit with parameters listed in Table ([T]) the first port is the qubit port. 
Here are the extracted Belevitch transformer matrices for each stage: 


Ti = 


T, = 


T., = 


T4 = 


T, = 


T« = 


T7 = 


/ -1.0000 

-0.0001 

-0.0010 

0.0008 

-0.7148 

-0.6993 

\ 0.0007 

0.6993 

-0.7148 

0.8933 

-0.0132 

-0.4493 \ 

-0.0132 

-0.9999 

0.0032 

-0.4493 

0.0030 

-0.8934 ) 

( 0.4315 

0.0060 - 

-0.9021 \ 

0.0127 

0.9998 

0.0127 

^ 0.9020 - 

-0.0169 

0.4314 j 

/ 0.0000 

-1.0000 

0.0030 \ 

1.0000 

0.0000 

0.0000 

\ 0.0000 

0.0030 

1.0000 j 

0.0000 

0.4254 

-0.9050 \ 

0.0403 

-0.9043 

-0.4250 

-0.9992 

-0.0365 

-0.0171 j 

-0.0416 

0.0024 

-0.9991 \ 

0.9299 

0.3659 

-0.0378 

0.3655 

-0.9306 

-0.0174 j 

-0.0006 

-0.9994 

-0.0342 \ 

1.0000 

-0.0006 

-0.0000 

0.0000 

-0.0342 

52 

0.9994 j 


(4.39) 

(4.40) 

(4.41) 

(4.42) 

(4.43) 

(4.44) 

(4.45) 








Figure 27: Fundamental mode (the TF'lOl mode) of the cavity with frequency /teioi = Q-STbGHz. Green color indicates 
electric field regions of higher magnitude compared to blue regions. 


j 

rj {ft) 

Lj (nH) 

Cj (nF) 

tj = l/rij 

^3,2 

^3,3 

1 * 

0.0923 

0 

1.1953 X 10-4 

0 

0 

0 

2 

0.0471 

7.1890 X 10^ 

2.4523 X 10-^ 

0.9478 

-0.0008 

-0.0259 

3 

0.0973 

2.7674 

7.7198 X 10-4 

0.0986 

-0.0002 

0.2050 

4 

0.1063 

2.7113 

7.8675 X 10-4 

0.0971 

0.0003 

-0.0020 

5 

0.2136 

3.0283 X 10^ 

1.7701 X 10-^ 

0.9915 

0.0037 

0.0018 

6 

20.7896 

2.7344 X 10^ 

1.0464 X 10-“ 

0.7657 

0.0002 

0.3708 

7 

21.4619 

2.7500 X 10^ 

1.0416 X 10-*^ 

0.7508 

0 

-0.0222 

8 

26.6330 

2.4557 X 10^ 

6.2335 X 10-y 

0.9959 

9.65 X lO-'^ 

-2.311 X 10-4 

9 

4.7957 

4.9851 X 10^ 

2.0961 X 10-^ 

0.8408 

0.0002 

0.0122 

10 

30.5600 

4.6115 X 10^ 

2.2697 X 10-^ 

0.8409 

0.0007 

-0.0623 

11* 

84.5207 

0 

2.4178 X 10-'^ 

0 

0 

0 

12* 

88.4419 

0 

2.2673 X 10-^ 

0 

0 

0 


Ri = 1.0837 X 10^ i ?2 = 1-1306 x 10^ R 3 = 7.7537 x 10^ 


Table 1: 3-port Brune circuit parameter values for the dataset corresponding to the setup in Figs. ii^ and (Eg. Stages 
marked with * are capactive degenerate stages. We note that in degenerate stages there are no n-type and i^-type transformer 
couplings. The first port is the qubit port. 
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T, 


8 — 

T9 = 


Tio = 


Til = 


Ti2 = 


Ti3 = 


0.9975 

0.0341 

-0.0615 

-0.0308 

0.9981 

0.0538 

0.0632 

-0.0517 

0.9967 

-0.9976 

0.0011 

-0.0685 

0.0032 

0.9995 

-0.0299 

0.0685 

-0.0301 

-0.9972 

-0.0011 

-0.9999 

-0.0109 

1.0000 

-0.0011 

0.0003 

-0.0003 

-0.0109 

0.9999 

-0.9775 

-0.1067 

-0.1820 

-0.1088 

0.9941 

0.0015 

0.1808 

0.0212 

-0.9833 

-0.0081 

0.9876 

-0.1566 

-1.0000 

-0.0080 

0.0013 

0 

0.1566 

0.9877 

-0.0978 

0.9951 

-0.0116 

0.9952 

0.0978 

-0.0001 

0.0011 

-0.0116 

-0.9999 


(4.46) 

(4.47) 

(4.48) 

(4.49) 

(4.50) 

(4.51) 


We note that all the Belevitch transformer matrices above are orthogonal - as expected - up to numerical 
noise. We would like to also mention that Tj being proportional to a permutation matrix for some j, 
1 < j < Ilf would re-order the ports without introducing a coupling between them. Note that all the Tj 
matrices, except for the first one, are roughly approximated by permutation matrices; this can taken to 
mean that there is only one “channel” whereby there is significant coupling among the three ports due to 
the first Belevitch transformer Ti. 



Figure 28: Multiport Brune stage circuit in the case of a non-reciprocal response. We observe the appearance of a multiport 
gyrator coupling the first port to the remaining ports with a gyration vector 7 as shown in H- 


54 





































5. Non-reciprocal Brune Stage 


The multiport Brune’s Algorithm described in Section ([3]) produces reciprocal stages for a reciprocal 
impedance response Z(s). If the response is non-reciprocal, i.e. if the blackbox in Fig. © contains 
circulators for example, the multiport Brune circuit extracted at each stage is slightly modified as 
shown in Fig. (1331) . [I^ . In Fig. (1331) we see that a multiport gyrator with a gyration vector 7 is extracted 
right after the resistance ri. The circuit symbol for this multiport gyrator is shown in Fig. (13^. It has a 
single port on the left and {N — 1) ports on the right with the following constitutive relations, |18l |: 


Vi 

V2 




(5.1) 


where Vi, Ii are the voltage and the current of the left port and V 2 , I 2 are vectors of length {N — 1) holding 
the voltages and currents of the ports on the right, respectively. 


h 


+ 

Fi 



7 



Figure 29: The multiport gyrator appearing in the non-reciprocal Brune stage in Fig. II28II . The gyration ratios are given by 
the vector y. 


The network Ni has the same time-evolution description as in Eq. (15.271) . However the input-output 
relation is slightly different with the appearance of the gyration vector 7 in the Di matrix: 


/ T)2AAlnj - 7 ^ -f T)2AB/ni 
7-|-D2BA/tll D 2 BB 


(5.2) 


That is 7 is extracted by taking the anti-symmetric part of Di. 

We note that the quantization of the multiport Brune circuit with gyrators - in the most general non¬ 
reciprocal case - is not solved by the formalism reported here, and remains an open problem. 


6. Conclusion 

As the scaling up of superconducting quantum processors continue, it is clear that one will need to deal 
with microwave circuits of increasing complexity. In fact, a route to scaleup is clearly taking shape: it is 
increasingly accepted that a clear step towards realizing a fault-tolerant quantum computing is by using 
the surface code . The surface code is a quantum error-correcting code with a relatively high threshold 
(estimated threshold error rate around 1 %) and with realistic requirements for, e.g., local couplings beween 
qubits. The implementation of the surface code will require maintaining the excellent fidelity of qubit oper¬ 
ations, made possible by high Q-factor resonators and other components, in systems of increasing scale and 
complexity |55l | . There will clearly be a need for highly accurate characterization of these multiqubit systems, 
so that their optimization can proceed according to a rational plan. Already in previous work, we have seen 
first attempts at such design efforts, with the focus on the problems with spurious modes that inevitably 
appear in complex microwave structures. In the significant contribution of off-resonant modes in the 
spontaneous emission rates of qubits was already demonstrated. This has led to a method for controlling 
the spontaneous emission rates of qubits by incorporating microwave filters in the superconducting quantum 
processors. There has been very recent research activity in this direction (H^. 1^: a new push is being made 
to design “Purcell filters” for the suppression of spontaneous emission of the qubits while keeping the qubit- 
resonator coupling at a strength required for measurement. Indeed, having multiple modes is not always 
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a nuisance. There are indications that one can engineer the multi-mode structure of microwave circuits to 
tailor qubit interactions via the help of filters for the realization of high-fidelity two-qubits gates . 

We are hopeful that the multiport method that we have provided here will be useful in the further 
efforts to achieve a more accurate analysis and design of these filter structures mentioned above, and of 
all other aspects of the structures needed to achieve surface-code action. Our work provides a complete 
algorithm for going from a proposed electric structure to its Hamiltonian description. It should be pointed 
out, however, that the present work leaves incomplete the one further step of this design process, in which 
the Schrodinger equation for this Hamiltonian is solved to obtain energy eigenvalues and eigenfunctions 
that are vital in, e.g, assessing the action of gating pulses, or of predicting qubit decoherence rates. Since 
our analysis endeavors to be exact, we obtain Hamiltonians with many degrees of freedom, with potential 
functions with a large hierarchy of stiffnesses . Such a high dimensional problem can be almost impossible 
to solve by straightforward discretization techniques. However, the multi-scale nature of the problem can 
be used to systematically identify “fast" and “slow" degrees of freedom in the problem [^; the principles of 
the Born-Oppenheimer approximation can be applied to have a very efficient approximate treatment of the 
fast variables, so that the dynamics of many fewer slow variables, which can typically be directly associated 
with the qubits or the principal resonant modes, can be determined accurately, leading to a physically 
appealing and computationally accurate description of the quantum physics of the system. We look forward 
to adapting the present mathematical analyses so that they are fully integrated into a full attack on the 
problem of scaling quantum hardware. 


7. Appendices 

7.1. Review of Lumped Element Circuit Quantization Methods 

In this appendix we will review two formalisms UM developed for the quantization of the lumped 
element superconducting circuits. We won’t attempt to do a full review of each formalism. We will rather 
content ourselves with describing how to combine parts of each formalism for the purpose of quantizing 
Brune circuits. In the following we will refer to 0 as “BKD” and to Q as “Burkard”; “KCL” stands for 
Kirchhoff’s current law and “KVL” stands for Kirchhoff’s voltage law. 

Both methods derive classical equations of motions for the lumped element circuits and identify canonical 
variables. They both start with a graph theoretical analysis by choosing a spanning tree in the circuit to 
write KCL and KVL relations involving current and voltage variables in the circuit. We will follow the graph 
analysis done in Burkard. This will allow us to write an equation of motion and to identify the canonical 
degrees of freedom of the circuit. We will need to make slight modifications to the theory to be able to treat 
Brune circuits. Once we have the equation of motion, we will interpret it as an equation of motion of BKD 
to do a dissipation analysis and compute relaxation rates. 

7.1.1. Derivation of the Equation of Motion by the Burkard’s Method 

The first step in Burkard is to find a spanning tree containing all the Josephson junctions, voltage sources 
and impedances in the circuit. One can also put linear inductors in the tree. However there should be no 
capacitors in the tree so that all capacitors are in the chord branches. Linear inductors are also allowed 
to be in the chord branches. Burkard assumes that there is no loop containing only Josephson junctions, 
voltage sources and impedances which is physically justified since in reality each loop would have a finite 
self-inductance. 

With the choice of such a spanning tree we can partition the current and voltage vectors as follows 


Itr = (Ij, li, ly, Iz) 

(7.1) 

Vtr = (Vj, Vl, Vy, Vz) 

(7.2) 

Ich = (Icj, T-c, lif) 

(7.3) 
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Vch = (Vc„ Vc, Vk) (7.4) 

where Itr and Ich are the vectors holding the currents through the tree and chord branches and Vtr and 
Vch are vectors holding the voltages across the tree and chord branches, respectively. Labels J, L, K, V, 
Z, Cj, C correspond to Josephson junctions, tree inductors, chord inductors, voltage sources, impedances, 
Josephson junction capacitances and ordinary capacitances, respectively. In terms of those loop variables 
Kirchhoff’s laws can be written as 


FIch = —Itr 


(7.5) 


F^Vtr = Vch - (7.6) 

where we have introduced the loop matrix F defined in Eq. (3) of Burkard. = (4>i,..., ‘hit) is the flux 
bias vector holding the external fluxes threading F fundamental loops of the circuit, each fundamental loop 
being defined by a chord branch. The loop matrix F can be partitioned as 


/ 1 

Fjc 

Fjk 

\ 

0 

Flc 

F lk 

0 

Fyc 

Fvk 


V 0 

F zc 

Fzk 

) 


(7.7) 


The first column is due to Josephson junction capacitances Cj’s shunting only the Josephson junctions. 

Burkard further assumes that the voltage sources and the impedances are not inductively shunted in the 
sense that 


Fvic = Fzk = 0 


(7.8) 


Then by writing KCL for Josephson junctions and tree inductors and KVL for chord capacitors Burkard 
derives the first-order equation of motion (Eq. (19) in Burkard - we fixed a sign typo) 


= Q - CyVy + J-cCz * Vc (7.9) 

where the vector holds the flux degrees of freedom corresponding to the fluxes across the Josephson 
junctions (J) and tree inductor branches (L) 


$ = 




(7.10) 


with €>,7 = $o<Pj/27r, tpj being the vector of phases across the Josephson junctions. The canonical charge 
variables are given by the vector Q 


Q = - 


Qj 

Ql 


— FkQk 


with 


-Ftc = 


Fjk 
F lk 


rc = 

and the capacitance matrices in Eq. (EH) are given by 

C = 


Fjc 
F LC 


Cj 

0 




(7.11) 


(7.12) 
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(7.13) 




Cv = J^cCFlc 


(7.14) 


Cz (w) = (^w) CFicZ (w) FzcC (7.15) 

where C is the diagonal matrix of ordinary capacitances in the circuit except Josephson junctions’ capaci¬ 
tances such that 

Qc = CVc (7.16) 

The last term J-cCz*^c in Eq. (|7.9D is the dissipative term. This term can also be written in frequency 
domain as 

J^cCz (w) Vc = - J-cCFicZ (w) Iz (7.17) 

\z can be written as function of flux degrees of freedom and voltage sources Vy as 

\z = -{iuj)[l + {iuj)FzcCFlc'^{oo)Y^FzcCT'E^ 

- (zee) [I + (zee) FzcCFicZ (ee)] FzcCF^c^y (7.18) 

If we define 

m = J^cCFIc 

rny = FvcCF^fz 

and as in Eq. (28) of Burkard (with a sign change) 

Cz (w) = - (zee) Z (ee) [l + (iee) FzcCFY'^ (ee)] (7.21) 

We can rewrite Eq. (|7.18|) in terms of the newly defined quantities in Eqs. (17.191) . (|7.20|) and (17.211) as 

Iz = Z^i (ee) Cz (ee) (ee) Cz (ee) (7.22) 

Replacing the solution for Iz in Eq. (17.221) in Eq. (17.171) we get 

^eCz (ee) Vc = —ihCz (ee) m^4> — ihCz (ee) rfiy Vy (7.23) 

Defining also 

Cz (ee) = ihCz (ee) (7.24) 

and 

Cv (ee) = rhCz (ee) rny (7.25) 

We can write the dissipative term in Eq. (17.231) in terms of the quantities in Eqs. (17.241) and (17.251) as 

J-cCz(ee)Vc = -Cz(ee)i-Cy(ee)Vy (7.26) 

Note that we have extracted an additional non-dissipative term Cy (ee) Vy in the equation above. 

Plugging the dissipative term in Eq. (17.261) back in the equation of motion in Eq. (17.91) we get in time 
domain 


(7.19) 

(7.20) 


(C + Cz (t)) * $ = Q - (Cy + Cy (t)) * Vy (t) 
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(7.27) 






















Note that the Eq. (17.271) is more general than Burkard’s equation of motioir in Eq. (25) of Q since 
we allow a general voltage source term Vy (t) with possible frequency components both at DC and AC. 
The vector Cy (t) is due to AC compoirents of Vy (t). In the case of a DC voltage source term Vy, 
Cy (w) Vy (w) = 0 for uj 0 and Cy (w) = 0 for a; = 0 so we recover Burkard’s equation of motioir. Eq. 
(17.2711 will lead to the more general Hamiltonian in Eq. (17.501) below. 

Taking the time derivative of both sides in Eq. (17.271) we obtain 


{C+Cz)*^ = Q-Cv*'Vv (7.28) 

Using the identity (C +Cz) * ^ = (^C +Cz'^ * $ it is interesting to note here that since the vector $ 
is equal to the vector of voltages across the Josephson junction branches and tree inductors, the factor 
(^C +Cz'^ is an admittance matrix and (^C +Cz'^ * $ = Q — Cy * Vy is the KCL for the tree branches. 
Taking the dissipative term C^ * $ in Eq. (17.281) to the right-hand side by noting the identity 

Cz*^ = Cz*^ (7.29) 

And using the Eq. (29) Q = — of Burkard with the potential U 

U ($) = -Lj^cos^J + -f- (7.30) 

where 


Mo = 


(7.31) 


N = (0 Ik f 

with being the inverse inductance matrix such that 



we get 

C * $ = — Icsirup — Mo$ — Cz * ^ — Cv * Vy — 


(7.32) 


(7.33) 


(7.34) 


(7.35) 


where U = is the diagonal matrix of critical currents. 

We would like mention here that the dissipation analysis that follows is to be applied only to unexcited 
systems, i.e. to systems where the voltage sources Vy (t) have only a DC component. 


7.1.2. Treatment of Resistors and Relaxation Rate Calculations Using the BKD Method 

Before writing a Hamiltonian for the equation of motion in Eq. (17.351) we have one more step to do. 
We again need to stretch Burkard theory slightly to include the resistors {i?i,..., i?Ar} shunting the last 
stage of the multiport Brune circuit in our analysis. The trick we used for that purpose as described in 
the main text for the one-port Brune circuit was to replace those resistors by capacitors and to make the 
substitution Cj ^ for 1 < j < V after we obtained the equation of motion. To do a similar analysis 
for the multiport Brune circuit, here we extend Burkard’s method to include chord resistors. 

We start by partitioning the matrix C 
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c = 


(7.36) 


{Co 0 \ 

V 0 Cfl y 

where Cq and C/{ are diagonal matrices of size M x M and N x N, respectively. Diagonal entries of Cq hold 
the values of ordinary capacitors in the circuit whereas is an auxiliary matrix which will be substituted 
later as 


(7.37) 

where R is the N x N diagonal matrix holding the values of the resistances shunting the last stage 

/ Ri 0 \ 

R= (7.38) 

\ 0 Rn / 

We also partition Tc respecting the partitioning of C in Eq. (17.3611 as follows 


= ( J^Co ^Cr ) 

where Rcqj ^Cr are submatrices with M and N columns, respectively. 

Now we can decompose the capacitance matrix C in Eq. (17.131) using Eqs. (|7.36ll and (17.391) as 


(7.39) 


where we defined 


C = 


Cj 0 
0 0 


Cj i 

0 I 

= Co + Cr 


TcCTI 


Tc,CoTI,\TcrCrTI^ 


(7.40) 

(7.41) 

(7.42) 


Co = 


Cj 0 
0 0 


Cr = TcrCrTc^ 


^Co Co^i 


Co 


(7.43) 

(7.44) 


Making the substitution in Eq. (17.371) we get the following dissipative term on left-hand side of the 
equation of motion in Eq. (17.351) 


Cr*^ = 

= (J-c„R-'i 

= 


(7.45) 

(7.46) 

(7.47) 


where in the first line we used differentiation property of the convolution operator, in the second line we 
used the Eqs. (|7.37l) and (I7.44|) and in the third line we defined 


= FcrR-^RI^ (7.48) 

which is frequency independent, justifying dropping the convolution operator in Eq. (17.471) . Now taking the 
dissipative term in Eq. (17.471) due to the last shunt resistors R to the right-hand side of the equation of 
motion in Eq. (17.351) and noting also the decomposition in Eq. (17.421) . we re-write the equation of motion 

(17351) 
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Co * $ = —Icsinyi — 7?. — Mo$ — Cz * ^ — Cv * Vy — N€>a; (7.49) 

Following Burkard we can write the following Hamiltonian for the lossless part of the dynamics described 
by the Eq. (j7.49l) 

77s = i (Q - (Cy + Cy) * Vy {t)f (Q - (Cy + Cy) * Vy (t)) + U ($) (7.50) 

where the potential function U {^) is defined in Eq. (17.301) . A point to note here is that although Cy 
defined in Eq. (17.141) is frequency independent, after the substitution in Eq. (17.371) it will acquire frequency 
dependence which we will analyze further down below in Section (17.1.31) . 

Comparing again the equation motion in Eq. (61) of BKD with the Eq. (17.491) above we identify 

^ (7-51) 

where Rbkd is now the resistivity matrix defined in Eq. (25) of BKD. We note however that is 

diagonal whereas TZ~^ is in general non-diagonal. 

To compute the contribution of the shunt resistor Rj alone to the relaxation rate we set Rk = oo for 
1 < fc < and k ^ j to write the term 77.“^$ in the equation of motion in Eq. (17.491) in the frequency 
domain as 


(7.52) 

where J^r^^Cr is the column of the matrix Tcr corresponding to the resistor Rj (or the capacitor Cr.). 
Comparing now 


Md (w) o (iuj)R^ ^■Frj,Cr-Fr.^Cr 

(7.53) 

where Md (w) is the dissipation matrix defined in Eq. (72)-(75) of BKD 

we identify the coupling vector 

mR^ = Tr.^Cr 

(7.54) 

and the dissipation kernel is given by 



(7.55) 


Using Eq. (93) of BKD we define the spectral density of the bath corresponding to the resistor Rj as 
(correcting the sign and dropping the scale factor) 


Jrj (w) = Im [Kr. (w)] 


UJ 



(7.56) 

(7.57) 


Hence the contribution of Rj to the relaxation rate is given by the formula in Eq. (124) of BKD as 

= l K0|fl.«, ■ *1 l>t*, (-o.)coth (7.58) 

where woi is the qubit frequency. Here we used flux variable and Hirj since we dropped the scale factor 
in Eq. (|7.56|) (BKD uses phase variable tp and the normalized vector m). 

We need to also consider the effect of last shunt resistors on the matrices Cz and Cy appearing in the 
equation of motion in Eq. (17.491) above. For this we first make the following partitioning for the matrix Fzc 


Fzc = ( FzCo FzCr ) 
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(7.59) 






















where the submatrices FzCo have M and N columns, respectively. 

Then using also the partitioning in Eq. (17.361) and the substitution Eq. (17.371) we can re-write the matrix 
Cz (w) defined in Eq. (I7.21|) as Cz,_r (w) 


CzM (w) 


where we have defined 


— (iuj) Zn (w) [l -I- {iuj)FzCoCoF^Co'^R (^)] 


(7.60) 

(7.61) 


with 


Zfi(c^) = Z[l + F2C„R-'Fic«Z] 


(7.62) 


Z = 


( ri 

V 0 



(7.63) 


being the diagonal matrix of in-series resistances ri,..., tm in the circuit. 

Since we will consider resistors one at a time we will take the limit R —> oo (which corresponds to open 
circuiting the last stage in the Brune circuit) in Eq. (17.621) above to define 


Cz..(c^) = (7.64) 

= -(*c^)Z[l + (*..)F2CoCoFicoZ]'' (7.65) 

We will see later below that C>z,r (w) is proportional to the dissipation kernel due to in-series resistors. 

We need to also update the coupling matrix rh defined in Eqs. (17.191) above to account for the effect of 
shunt resistors on the terms due to in-series resistors. First we again make a partitioning 

Tc = ( Tc, Tcn ) (7.66) 

Then using this partitioning and Eqs. (17.191) . (I7.36L (17.371) and (17.591) we can write the decomposition 


where we defined 

m = mo -1- m_R (w) 

(7.67) 


“lo = .F Co^oF^^Cq 

(7.68) 

Here, however. 

iTifl (w) = J'cbR’^fIc'b 

we will onlv use the zeroth order term in Ea. (17.681) to define 

(7.69) 


Cz,r (w) = moC^^r (w) HiQ 

(7.70) 


Note that the frequency dependent factors (iuj) ^ appearing in higher orders terms due to the coupling 
matrix (w) in Eq. (17.691) above can be absorbed in the corresponding Cz (w) matrices if one wants to 
investigate higher order effects. 

Comparing Eq. (17.491) with the equation of motion of BKD (Eq. (61) in [^) we make the following 
identification 


Md o Cz,r 
Mci (w) O —U}^Cz,r (w) 
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(7.71) 

(7.72) 



























where on the left we have the dissipation matrix defined in BKD and on the right we have the dissipation 
matrix appearing in Eq. (17.491) . The capacitance matrix Cq in Eq. (17.491) maps directly to the capacitance 
matrix C in BKD and the stiffness matrices Mg appearing in Eq. (I7.49|) and BKD map to each other. 

To do a dissipation analysis we will treat the Eq. (17.491) as an equation of motion of BKD and we will 
do a Caldeira-Leggett analysis using the identification in Eq. (17.7211 . By Eq. (64) in BKD we have 


{uj) = (w) 

And bv Eq. (17.701) 

(7.73) 

Cz,r (w) = (w) 

Hence bv usine Ea. (17.721) we further identify 

(7.74) 

(^) (^) 

Coupling matrices m and mo map to each other directly 

(7.75) 

rn O trio 

(7.76) 

To compute the contribution of each resistor to the relaxation rate l/Ti 
seperately. In that case by Eq. (73) in BKD we have 

we will treat each resistor 


K (w) = L/ (w) 

where K (lj) is a scalar. Hence bv Ea. (17.751) we eet 

(7.77) 

/O (w) = 

for the eauation of motion in Eq. (17.491). where we 

— Lo'^Cz,rj (w) 

defined 

(7.78) 

Czr.i(^) — ^Zr(^^\ n £ » / - 

^’^3 y > \ yirk=0,fork^j 

(7.79) 

= —iUJTj 

1 -1- ILUrjFrj^Co^oFr^^Co 

(7.80) 


where Fr^^Co i® ^^le row of the matrix FzCo defined in Eq. (17.591) . We note that Cz,rj (w) is a scalar. 
That is we are treating only the in-series resistor rj by short circuiting the other in-series resistors setting 
rfc = 0 for fe j. Note that shunt resistors R are already open circuited by taking the limit of R —?► c» in the 
definition of Cz,r (w) in Eq. (I7.64|l . We also define the coupling vector m.j of the bath due to the resistor 
Tj to the circuit degrees of freedom by taking the column of the coupling matrix rng in Eq. (17.6811 . 

Using Eq. (93) of BKD we define the spectral density of the bath corresponding to the resistor rj as 
(again correcting the sign and dropping the scale factor) 

Jj (w) = Im [Kj (w)] (7.81) 

We can now write the contribution of the resistor rj to the relaxation rate using the formula in Eq. (124) 
of BKD as 


^ 1(0 |mj • $1 1)1^ Jj (wQi) coth (7-82) 

where wgi is the qubit frequency. Here we again used flux variable $ and rhj since we dropped the scale 
factor in Eq. (17.811) . 
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7.1.3. Analysis of Voltage Source Couplings 

Now we will update the Cy and the Cy matrices to account for the effect of the resistors R shunting the 
last Brune stage. Making first the partitioning 

Fyc = ( FyCo FyCfi ) (7.83) 

Then using the definition for Cy in Eq. (|7.14|) together with Eqs. (|7.36l) . (I7.37|l . (17.391) and (17.831) we can 
decompose Cy as 


Cy = Cy.o + CvM (w) (7.84) 

where we have defined 

Cv,o = J^CoCoF^Co (7-85) 

and 

Cy.K (w) = (*w)-' J^Cn'R^-^F^Cn (7-86) 

We note the frequency dependence of Cv,r. 

A similar analysis can be done also for the rhy matrix defined in Eq. (17.201) . Then again using the 
partitioning in Eq. (|7.83l) and Eqs. (17.201) . (17.361) . (|7.37l) and (17.591) we can write the decomposition 


my = myo + my_R (w) (7.87) 

where we defined 

mv.o = FycoCgF^Co (7.88) 

mv.R (w) = (fw)"^ Fyc„R~^F|c'^ (7.89) 

Now using the definition for Cy in Eq. (17.251) together with the decompositions Eq. (17.671) . (|7.87l) and 
the updated Cz (w) defined as Cz,r (w) in Eq. (|7.60l) we can re-write Cy as 


Cy(w) = (7.90) 

= (riio-f iTifl (w)) (w) (rhyo 4-my/i (w))^ (7.91) 

Unlike the dissipation matrix Cz we will keep the full expression for the term Cy (w) since it is not a 
dissipative term. However we will combine the frequency dependent term Cy n (cu) defined in Eq. (I7.86|) 
with Cy (w) in Eq. (17.911) to define 


Cy,R (oj) = CyR (uj) + Cy (uj) 


(7.92) 


Using the matrices defined in Eq. (17.701) and (17.921) above we write again the equation of motion in Eq. 
(I7.49|) as 


Co * ^ = —Icsinyi — 7Z — Mo$ — Cz,r * ^ — Cy,R * Vy — (7.93) 

The Hamiltonian of the system described by the equation of motion in Eq. (17.931) is given by 

Hs = i [Q - (Cyo + CyR (t)) * Vy (t)f C^^ [Q - (Cyo + CyR (t)) * Vy (t)] + U ($) (7.94) 

where the potential function U ($) is defined in Eq. (17.301) . This time-dependent Hamiltonian is the 
extension of the Hamiltonian given in Eq. (36) of Burkard to a time-dependent voltage source vector 
Vy it). 
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1.2. Derivation of the Effective Kirchhojf’s Voltage Law 

7.2.1. Effective Kirchhoff’s Voltage Law for the One-Port State-Space Brune Cireuit 

Here we present the effective Kirchhoff analysis for the Kirchhoff’s voltage law for the one-port state- 
space Brune circuit in Fig. ©• The treatment here will be along similar lines with the analysis we did in 
Section (1^ to get the effective Kirchhoff’s current law. 

We start by writing the Kirchhoff’s voltage law for the circuit in Fig. 0 


F^VtF = Veh 


where the voltages are partitioned as 


Vtr= (K7,Vi,Vz,Vf^) 

(7.95) 

Vch= 

(7.96) 

and the loon matrix matrix F"^ is oartitioned bv Ea. (12.7011 as 


/ \ 

■^T _ 1 ^ JC ^ LC ^ ZC ^TC \ 

^ ~ \ WT -tT t^T 

\ ^ JT ^ LT ^ ZT ^ TT J 

We will show that we can get an effective Kirchhoff’s voltage law 



(FT)e//ve//^v:" 

for some effective loop matrix partitioned as 

(fT)'" = ( (FlaT" {Flaf> ) 

by eliminating the transformer branches’ voltage variables 


= {Vj,Vl,Vz) (7.97) 

KU = Vc (7.98) 

First we note the following in the circuit in Fig. © 

Hence by the transformer voltage relations we have 

(7.100) 

= - NVz, (7.101) 

where N is the diagonal turns ratio matrix defined in Eq. (I2.81|l . 

Now writing voltages of the chord capacitors as a function of voltages of tree branches in the circuit in 
Fig. © we get 


with 


Vc =F 


%Vj 


fIcVl 


■FzcVz + F, 


TC^T 


(tr) 


(7.102) 
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(7.103) 


■pT 

^ JC 


Fic 


zc 


IpT 

^TC 


/ 1 

V 1 

(1 

1 1 0 

1 1 ••• 
Vi 1 ••• 
(1 

1 1 0 


1 1 
Vi 1 
/ 0 
1 ■■ 

Vi •• 


1 / 


1 / 


0 \ 


0 
1 / 


(7.104) 


(7.105) 


(7.106) 


where is a vector of length (M + 1) and FT^, FX^, F!^^ are matrices of size (M + 1) x M. 
Using Eqs. (rmm and we get 

Vc = FjcK/ + {Flc - F?cN)Vi + F^^Vz 


(7.107) 


from which we conclude 


Uc)-" 

II 

(7.108) 

-pT \e// 
^ LC) 

= fIc-F?cN 

(7.109) 

f5c)'" 

II 

O 

(7.110) 


and the effective Kirchhoff’s voltage law 

(7-111) 

Comparing above Eqs. (I7.108l) - (|7.110l) with (|2.86|) . (12.881) and (I2.91|) of Section (12.61) we conclude 


[F^jcY" = 

(f'/j; 

T 

) 

(7.112) 

II 

(f7i: 

T 

) 

(7.113) 

fIcY" = 

(f7i; 

T 

) 

(7.114) 


Hence 

(jTyff ^ ^-peff^T (7.115) 
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1.2.2. Effective Kirchhoff’s Voltage Law for the Multiport Brune Circuit 

In this section we will derive an effective Kirchhoff’s voltage law for multiport Brnne circuit in Fig. m- 
We write Kirchhoff’s voltage law 


= V, 


ch 


(7.116) 


with partioning in Eqs. (13.491) and (|3.50|) for voltage vectors Vtr and Vch respectively. We note that the 
relation in Eq. (|7.116|1 derives from a graph theoretical analysis Q of the multiport Brune network. 


where 


and 


vectors 

and V 

for the transformer branches 

as follows 

\ T — 



(7.117) 

< 

II 



(7.118) 

— 

^ n 


■ • 5 ^ nu j 

(7.119) 

V T — 



(7.120) 

— 



(7.121) 

Y^^') — 

^ n 

{ykW- 


(7.122) 

V rp — 



(7.123) 

yW — 

* V 


•■><) 

(7.124) 


with 


V, 


{L){R) 


v(«) = 


/ \ 

V j 

( \ 


\ / 


(7.125) 


(7.126) 


where are vectors of length for 1 < j < M + 1 and are vectors of length (N — 1) for 

1 < j < M- 

We will move in the opposite direction to the direction we have chosen in Section (13.21) through the 
multiport Brune circuit while deriving effective loop matrices. That is we will start at the leftmost part of 
the multiport Brune circuit in Fig. m and move to the right. We again assume that shunt resistors Rj’s 
in the last stage are replaced by capacitors Cr. ’s for 1 < j < as shown in Fig. (EOl) and that all of the 
ports of the multiport Brune circuit are shunted by Josephson junctions as shown in Fig. m-, we get the 
relation 
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(7.127) 

where Vj is the vector holding the voltages across the Josephson junctions shunting the ports of the multiport 
Brune circuit. We again assume the following ordering for the capacitors 


(7.128) 

The voltages of inter-stage transformers are given by Eq. (18.6911 

(7.129) 

for 1 < j < M -I- 1, where Tj is the {N x N) Belevitch transformer matrix of the stage. The voltages of 
consecutive inter-stage transformers are related by 


= AjeiW, + vJV,^ + AjV^f 


(7.130) 


T 

for 1 < j < M, where ei = ( 1 0 ... 0 ) is the unit vector, Aj is the {N x N) matrix defined in Eq. 

(13.7811 and Vj is the row vector defined in Eq. (13.8011 . 

We can write the voltage Vc^ across the capacitor Cj at the stage as 


Vc,=VL,+Vr,+e^V. 


(fl) 


(7.131) 


We are now going to iterate over the index j starting at j = 1 through the stages in the multiport Brune 
circuit using the Eqs. (I7.127L (17.12911 . (17.13011 and (17.13111 


V, 


(H) 


V, 


(L) 


T2 


V, 


(R) 


V, 


{L) 


yc2 


= TfVj 

AfeiW, = AfeiWi TAfTfVj 

Vl, + = Vl, + W, -b efTfVj 

=T^AfeiW, TT^AfTfVj 

AlTlAle^Vr,+Ale^Vr,+AlTlvlVL,+vlVL, + 
TA^T^AfTfVj 

elTlAle^Vr,+Vr,+BlTlvlVL,+VL, + 

+e^T^AfTfVj 
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^C3 


Tjvg 

T^A^T^AfeiK-, + 

+T J A^ T^ v'IVl,+ T^ vl Vl, + T^ A^ T^ Af Tf V j 
Ai’eiK3+t.Jyz.3+Ai’vg) 


A^ T^ A^ T^ Af eiV;, + A^ T^ A^ eiK, + A^ 01^3 + 

+Al T J Al Tl vIVl,+ a J vl Vl, + yL3 + T J A^ T^ Af Tf V j 


Vz. 3 +K 3 +ervg^ 

efTjA^TjAfeiK^ +e?^TjA^eiK3 +K3 + 

+ef Ti’ A^ Tl vf Vl, + ef Vl, + Vl, + 

AefTi’A^T^AfTfVj 


Hence from the above relations we define the following for 1 < j < M and 1 < k < M 


(FJc)'" U,k) = 0 


{Flaf^ij,k) = ejTl^,vl 

f- _ ^TrrT xT rpT 




{j,k) = 0 

{F^cY" ij^k) = l 

, (j, fc) = ef TjAj_,Tjli... 


for j < k < M 
for k = j 
for k = j — 1 
for 1 < fc < j — 2 

for j < k < M 
for k=j 
A^ei for I < k < j — 1 


and for 1 < fc < A 


I (J> k) = efTf efc for j = 1 

\ {F^jcY^^ (J: k) = ejTjAj_,Tj_, ... AjTfe, forl<j<M 
To compute effective loop submatrices for j > M we note from Eq. (17.1321) the following 


= AlTl...AiT^Afe,Vr, + ... + Al,e,Vr^ + 
+Al,Tl ...A^T^v{Vl, + ... + vIVl^ + 
+AlTl...AjTlYj 

Writing also the voltage relations for the last stage transformer using Eq. (|7.129l) 


And noting 


V 


(R) 

Tm+1 


= T 


T 

M+1 


V 


(L) 

T^f+i 


We conclude using Eq. (17.13611 


Vch 


V 


(«) 

Tm+1 


rpT 

^M+1 


V 


{L) 

Tm+1 
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(7.132) 

(7.133) 

(7.134) 

(7.135) 

(7.136) 

(7.137) 

(7.138) 

(7.139) 










(7.140) 


Vch — T^_|_;^A^T^ ... A^T2 A^eil4i ... + T^_,_jA^eil4M + 
+T^+i A^T^ ...AItWiVl,... + + 

+Tl^,Al,Tl,...AlTlVj 

From Eq. (17.1401) above we define for 1 < fc < M 


(^Lc) 0> 

{^LcY^^ 0'^ = eJ_MTM+i^^M 


for k < M and M + 1 < j < M + N 
for k = M and M + l<j<M + N 


(7.141) 


{ {fIcT" 0) k) = eJ_^Tl^,AlTl ... Al^,Tl^,Ale, forM + l<j<M + N (7.142) 

where Bj is the unit vector of length N non-zero only at its entry such that (k) = 0 for fc ^ j and 

ej (j) = 1- 

And for 1 < fc < A 

{{F^jcT^^ ij,k) = eJ_j^Tl^,AlTl...AlTlek for PI + 1 < j < M + N (7.143) 

Hence from Eqs. (I7.133L (I7.134E (|7.135l) and (I7.141|l . (I7.142|) . (I7.143|l we conclude 

^pTyff ye// ^ yeff (7.144) 

with 

(rq'" = ( (Fje)'" (Fla)’" (Flo)’" ) I’-liS) 

and comparing Eqs. (I7.133|l . (I7.134|) . (17.13511 and (17.14111 . (I7.142E (I7.143|l to Eqs. (13.8311 . (13.8811 and (|3.91|1 

of Section (13.21) we conclude 

(fT)^// ^ (7.146) 
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